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We study the magnetic phases of two coupled two-dimensional electron gases in order to de- 
termine under what circumstances these phases may occur in real semiconductor quantum wells 
and what the experimental properties of the broken-symmetry ground states may be. Within the 
local-density-approximation to time-dependent density functional theory (DFT), we find a phase 
transition signaled by the vanishing of the intersubband spin-density excitations at low but accessi- 
ble (~ 10^" — 10^^ cm"'^) electron densities. Through a self-consistent Hartree-Fock calculation, we 
associate this transition with an antiferromagnetic phase and study the phase diagram, thermody- 
namics, and collective modes in it. The collective modes are in principle observable in inelastic light 
scattering experiments, and we discuss the implications of our calculations for these measurements. 
We also examine the ferromagnetic transition in both single and double quantum wells within the 
local-spin-density approximation to DFT and obtain a critical density which depends on the well 
width and which is far below that of the antiferromagnetic transition. 
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The study of spin instabilities holds a fascination for both theoretical and experimental condensed matter physicists. 



< 



. I. INTRODUCTION 

0^ 

I To theorists, these instabilities illustrate the qualitatively new states of matter which may result in simple systems 
. through the presence of electron-electron interactions. To experimentalists, they lead to novel phases with unique and 
\^ ' potentially useful properties. Consequently, the search for systems which exhibit new or unusual spin instabilities is 
On > an area of active research, and any guidance that theory can provide which suggests avenues for this investigation 
should be welcome. 

, One class of systems ready for serious exploration are single- or double-quantum- well structures at low density. The- 
oretically, the low-dimensionality of these structures restricts the phase space available for electron-electron scattering, 
increasing the relative importance of the interaction and thereby enhancing the potential for novel phase transitions. 
Another, more practical, reason for looking at these structures is that they can be physically realized in semiconduc- 
Q tor space-charge layers and, in particular, ultra-pure modulation-doped GaAs/Alj^Gai-^jAs heterostructures. These 
devices may be fabricated with precisely controlled dimensions that are nearly free of defects and which are tunable 
over a wide range of densities and band structures. This freedom yields a large parameter space in which interesting 
effects may be found and explored with both experimental and theoretical tools. 

In isolating the interesting regions of this parameter space, researchers can be guided by simple energetic considera- 



^. . . ^ 

jH . tions. When the kinetic energy of an electron gas dominates the Coulomb repulsion, as at high electron densities, the 
- - ■ system behaves like a nearly ideal Fermi gas. Therefore, novel phases will be found only in the regime in which the 
kinetic energy is smaller than or of the same order as the potential energy. Two ways of reaching this limit suggest 
themselves: applying a magnetic field or reducing the electron density. In a strong magnetic field, the kinetic energy 
is quenched by the Landau quantization, leading tp_a variety of strongly correlated quantum phases, the best-known 
example being the fractional quantum Hall liquid.EJcl At extremely low densities, the electrons prefer to crystallize, as 
predicted some time ago by Wigner.Q Experimental evidence for this crystallization has been somewhat equivocal in 
semiconductor heterostructures, but it is clearly seen in a two-dimensional electron gas suspended above the surfacej 
of liquid helium.E3 For somewhat higher densities, magnetic instabilities to spin-density wavetil or ferromagneticE3 
phases have been proposed. 

The parameter space defined by the energetics of a single low-dimensional electron gas is well defined and thoroughly 
explored. To find new effects, many current investigations of these systems add an additional degree of freedom by 
coupling two two-dimensional electron gases together. The resulting two-layer system has an additional energy scale 
due to the splitting of the isolated quantum well levels into symmetric and antisymmetric components which competes 
with the intra- and inter-layer Coulomb energies. In a high magnetic field, this additional degree of frefidgm can lead to 
the disappearance of odd-integer and the appearance-Qf even-integer fractional quantized Hall stepsl'EHZl and possibly 
to spontaneous charge transfer between the layers.Eja 
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In zero-field, a number of experimental results and theoretical predictions have also appeared in the literature 
Many of these zero-field studies focus on the behavior of the collective excitations involving transitions between the 
symmetric and antisymmetric levels and which are therefore unique to two-layer systems. These excitations can be 
produced in either the charge or the spin channels and have energies which are sensitive to the many-body interactions 
in the system. For example, in quantum wells with cither square or parabolic confining potentials, theoretical work 
suggests that there exists a critical density below which the many-body corrections to the charge-density excitations 
cause the energy of this mode to drop below the energy O^nthe symmetric-antisymmetric splitting.Ea This effect was 
recently observed by inelastic light scattering experiments.tj Another experimental study of coupled double-quantum- 
well systems reveals that the spin-density excitation energy abruptly merges with the continuum of intersubband 
single-particle excitations as the electron density is increased beyond the poiHt|at which the second-lowest subband 
begins to populate.Ell Time-dependent local-density approximation calculationsE3 agree quite well with this experiment 
for most of the range of densities used, but the abrupt merge seems to demand a more refined calculation. 

Of most relevance to the subject of this paper, recent calculations of the intersubband spin-density excitations in 
coupled double-quantum- well systems indicate that the energy of the lowest intersubband spin-density excitation may 
vanish at sufficiently low densityEJ The obvious interpretation of this collapse is that it indicates an electronic phase 
transition from the metallic Fermi-liquid phase to a condensate of zero-energy spin-density excitations. Since these 
excitations involve an intersubband electronic transition accompanied by a spin-flip, this condensate has been termed 
a spin-triplet intersubband exciton liquid. The word "exciton" here does not refer to the usual bound state of an 
electron in the conduction band and a hole in the valence band of the semiconductor; rather, it is used as a remindpc 
that the final-state interaction or vertex correction is included in the calculation of the spin-density excitations .l3 
In general terms, this transition involves the electronic spin in a fundamental way, suggesting that the new ground 
state would have non-trivial magnetic properties. Hence, the region of parameter space including the quantum well 
structures exhibiting this spin-density excitation collapse are ideal candidates for the study of novel spin instabilities. 

In this paper, we examine the question of spin instabilities in such single- and double-quantum-well structures in 
the absence of an external magnetic field. Our goal is to determine in what structures and under what conditions 
these instabilities may occur and what the experimental signatures of the-iiew ground states may be. Our primary 
interest is in the spin-density-excitation-collapsed phase discussed above,EJ but we also explore the more general 
question of ferromagnetism in single- and double-quantum-well structures. We predict that the SDE-softened phase 
will occur in fairly typical coupled double-quantum-well structures at low but accessible densities of order 10^*^ — 10^^ 
cm~^, and we show that this phase corresponds to antiferromagnetic order of the spin densities in the two wells. By 
constructing a minimal model of the antiferromagnetic state and treating it within mean-field theory, we are able to 
discuss the qualitative features of this phase. In particular, we find that the transition to this state may occur at 
temperatures on the order of the splitting between the lowest two subbands in the quantum wells, which can be around 
10 K. Moreover, although transport measurements will likely show no pronounced anomaly at the transition, both the 
collective excitations and the specific heat show distinctive features which can be used to identify the antiferromagnetic 
phase. Lowering the electron density further, our calculations indicate that these systems re-enter the normal state 
and then enter a ferromagnetic phase at densities around 10^ cm~^. In wide single-quantum- well structures as well, 
both ferro- and antiferromagnetic phases exist, and we study the critical density of the ferromagnetic transition as a 
function of the well width. As expected from the increasing importance of exchange effects in lower dimensions, we 
find that the ferromagnetic phase is stable below a critical density which increases as the well width decreases. 

This paper is structured in the following way. In Sec. H, we discuss the formalism which underlies our computations. 
In the normal and ferromagnetic states, this formalism is density functional theory in the local- (LDA) and local- 



spin-density approximations (LSDA), respectively. A time-dependent version of the LDA is also reviewed in Sec. [I A , 
which provides quan titative results for the collective excitation spectra that we discuss in Sec. |III^ Additionally, we 
describe in Sec. II B the equations for the self-energy and the density response function in a self-consistent Hartree- 



Fock theory, which is the basis of our calculations within the antiferromagnetic state. In Sec. Ill, we reproduce 



and extend the results of Ref. which predicted the softening of the intersubband s pin -density excitations, and 
characterize the phase transition on the paramagnetic side in greater detail. Section [V contains a study of the 
ferroma gne tic transition in these systems within the LSDA, and we demonstrate that that the instability predicted 



Sec. 



Ill occurs at much higher density than the ferromagnetic transition and so cannot be associated with it. 
We also take a first step in the study of the ground-state spin-polarization of the inhomogeneous electron gas in 
semiconductor quantum wells by calculating the critical density of the ferromagnetic transition in square single 
quantum wells as a function of well width at zero temperature. In Sec. we develop a simple model for the spin- 



density instability of Sec. Ill and treat it within the self-consistent Hartree-Fock theory of Sec. [IB. We study the 



ground state, thermodynamic quantities, and the collective mode spectrum in the broken-symmetry phase and discuss 



their experimental ramifications. Section VI presents some speculations on the importance of non-trivial spin-density 



modulations transverse to the layering direction in these quantum well structures and summarizes the results of this 
paper. 
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II. FORMALISM 



A. Density-Functional Theory 

1. Unpolarized Electron Gas (LDA) 
The central aim of this Subsection is to develop the formalism used to compute the intersubba nd c ollective exci- 



tations in the unpolarized state of double and wide single quantum wells that we will use in Sec. III. We compute 
these collective excitations within the so-called time-dependent local-density approximation (TDLDA) . This approach 
was first employed fftijAndoEII to compute intersubband charge-density excitations (CDEs) and later extended by 
Katayama and AndcE3 to study resonant inelastic light scattering in semiconductor structures. The use of inelastic 
light scattering in these systems is motivated by the fact the charge- and the spin-density excitations (SDEs) cooinla 
to the light polarization differently, and this fact allows a selective measurement of both types of collective modesE3't£l 
Detailed desceiBtions of the TDLDA method for calculating CDE and SDK energies and spectra have been given in 

we 



the literatureBaH However, for the sake of completeness and to facilitate the discussion of our results in Sec. Ill, 
describe the TDLDA approach in some detail in the remainder of this Subsection. In the next Subsection, we shall 
generalize this formalism to allow for the possibility of a spin-polarized ground state of the electron gas. 

The first step in the TDLDA calculation of the intersubband excitations consists of obtai^ng,the renormalized 
subband energies in the local-density approximation (LDA) of Hohenberg, Kohn, and Sham.EirL3 We choose a co- 
ordinate system with the z-axis is along the direction of confinement in the quantum well structure. The effective 
single-particle Schrodinger equation or Kohn-Sham equation for this system then reads 



-— V^ + VbffW *(R) = £;*(R), (1) 
Im* J 

where we have assumed that the effective electron mass m* is constant across the well, R = (r, z) denotes a three- 
dimensional vector, and the self-consistent effective potential Veff{z) is given below. The in-plane {xy) and z- 
dependences can be separated, and, due to the assumed translational invariance in the xy plane (the localized donor 
charges are assumed to be smeared out uniformly in the plane), the cigcnenergies and wave functions become 



2m* 



vE'„k(R) = 4Te""""'^"W- (3) 
V A 



and 



In these expressions, A is the sample area, k is the in-plane wave vector of the electron, and e„ and 0„(z) are the 
solutions to the one-dimensional Kohn-Sham equation 



+ VeFf{z) (t>n{z) = Sn(t>n{z). (4) 



2m* dz"^ 

The effective single-particle potential 

Veff{z) = Vconf{z) + Vh{z) + Vxc{z) (5) 

contains the confining potential of the bare quantum well Vconf{z), and the self-consistent Hartree and exchange- 
correlation potentials Vh{z) and Vxc{z), respectively. The Hartree potential takes into account the average electro- 
static interaction with the other electrons and the positively charged donor ions, and is given by the solution of the 
Poisson equation 

d'^VH(z) 47re^ . , , , ^ 

fj ' = —Hz)~Nd{z)], (6) 

where e is the static dielectric constant of GaAs, n{z) is the electron density, and Nf){z) is the density of positive 
donor charges, which are assumed to be located far from the quantum wells. Integrating Eq. (^) twice, one obtains 
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Vh{z)^ —i^J dz'{z- z')n{z') + z j dz'n{z')- ^zj, (7) 

where Ns is the eleptronic sheet density. For t he exc hange-correlation potential, we use the parametrization due to 
Ceperley and Aldeica given by Eq. (H) in Sec. |lIA2| for both the spin-polarized and -unpolarized cases. 



Once the subband energies e„ and wave functions (j)n{z) are obtained by solving Eq. (jj) numerically, the z-dependent 
electron density is calculated from 

n(z)=5,^/(i?„k) l*nk(R)|', (8) 

nk 

where the factor gs accounts for the spin degeneracy and f{E) = 1/ {e^^ -t- 1) is the Fermi-Dirac distribution function 
with (3 = \/T the inverse temperature {ks = 1 throughout this paper). This density may be rewritten as 

n(z) = ^n„ |(/)„(z)|^ (9) 

n 

with the subband occupancy rtn given by 

— 

A 



= (10) 



k 

The chemical potential is determined implicitly by the relation 



,11) 



At zero temperature, f{E) Q{—E) and n„ becomes 

Tin ^ gsNa{£F - £n)Q{£F - e-n), (12) 

where Q is the step function, £p \s the Fermi energy, and Nq = m* I^-kT? is the two-dimensional, single-spin density 
of states. 

The LDA electronic structure for the spin-unpolarized case is thus obtained by solving Eq. (^) together with Eq. (|^) 
for Vf/[n(z)], Eq. (^ ) for Vxc[?^(^)]i and Eq. (g) for n{z) self-consistently. The results of applying this procedure to 
a typical double-quantum-well structure are shown in Fig. |l|. 

From the wave functions and eigenenergies of the LDA calculation, we can compute the collective modes of the 
confined electron gas which are visible in inelastic light scattering experiments. The relation between the ccess section 
of resonant inelastic light scattering by electronic excitations was obtained by Hamilton and, McWhorterta for bulk 
systems was adapted to the case of semiconductor heterostructures by Katayama and Ando.Ea These authors showed 
that the cross section for inelastic light scattering by CDEs is proportional to the imaginary part of the reducible 
electronic polarizability function 11 with the proportionality factors depending on the details of the band structure 
of the host material. These factors vanish for perpendicular polarizations of the incoming and scattered light, and 
are maximized for parallel polarizations. Thus, CDE spectra are measured in practice within the so-called polarized 
configuration, i.e., with parallel polarizations of the two beams. On the other hand, the scattering cross section due 
to spin-density excitations is proportional to the imaginary part of the irreducible electronic polarizability function 
n, which is also called the spin-polarizability function, and contains prefactors that are maximized for perpendicular 
polarizations of the incident and scattered beams. This is the usual geometry employed in measurements of SDE 
spectra, and is referred to as the depolarized configuration. Since we are interested in the properties of the electron 
gas confined in the semiconductor structure, it is sufficient to calculate the electronic response functions 11 and 11 and 
ignore the band-structure-dependent factors. 

We compute these response functions within the TDLDA, which is equivalent to calculating the irreduciijle polar- 
izability function including a static, q-independent vertex correction in the ladder diagram approximation.Ej Within 
this approximation, the integral equation for the irreducible polarizability Il{q,Lu) can be solved exactly and gives 

n°(g,a;) 

"^'^'")^i + c/.cno(,,.)' ^''^ 

where Il^(q,uj) is the leading-order polarizability function and Uxc is the static, g-independent vertex function. From 
Dyson's equation for the effective Coulomb interaction,E3 one obtains the reducible polarizability function 
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nfa^)- , y (14) 

where Uh{(1) is the Fourier transform of the bare Coulomb interaction. Combining Eqs. ( [l3| ) and (14), we obtain 

m ^ n"(g,a;) 

"^'^'"^^i-(c/H(.)-t/xc)no(,,.)- (1^) 

In a confined electron gas system, where the confinement discretizes the single-earticle energy levels, the collective 
excitations must be calculated within a generalized dielectric function formalism.EZI In this context, the functions Uxc 
and Uh{(i) are replaced by matrices with indices labeling the different subbands. Within TDLDA, we have 

jxc _ [ [ .lJMJ.^MJ.^^^xc , 

and 



Uijl^nn^- dz dz'(^,{z)(^,{z)-^{z)5{z~z')cl)^[z')c^^{z') (16) 



U"^n.n{q) = ^/ dzdz'c^,{z^,{z)e-^\^-^'\ (bm{z')U^'), (17) 

where e is the background dielectric constant. 

The reducible polarizability function n(q, qz,uj), whose imaginary part is proportional to the spectrum of the CDEs 
and to the Raman intensity in the polarized configuration, is given byES 

n(q,(7.,a;) = J dz J dz'e-''^^^'~'"> Il{z, z';<i,oj) (18) 

with 

n(z,z';q,w)= ^ 0i(z)0j(z)ny^fc;(q,tj)0fe(z')(/);(z'), (19) 

i,j,k,l 

nij,fci(qi^) = n°(q,w)(5ifc(5j7 +^n^j(q,w)J7ij,m„((?)n™„,fci(q,t^), (20) 

m.n 

C/.,,™„(q) = t/.?™„(q) - U^l^„, (21) 

and 

n.(q> -) - 2 1. ^^.(k + q)) - E.ik)) - n{u + ■ ^''^ 

In these equations, subscripts are the subband indices, q and k are two-dimensional in-plane wave vectors, and (pj and 
£j are the LDA-calculated subband wave functions and energies. In addition, 11°^- is the leading-order polarizability 

function for the transition i — > j, Ej{k) — ej + f{E) is the Fermi factor, and 7 is a phenomenological inverse 

scattering time; at T = 0, an analytic expression for II"^ can be found in Ref. We note that the random-phase 
approximation (RPA) is obtained in the subband representation by removing the vertex correction U^^^ in Eq. (pl|). 

The imaginary part of the irreducible polarizability function 11 is proportional to the SDE spectrum and to the 
Raman-scattering intensity in the depolarized configuration. In the subband representation, the calculation of 11 is 
analogous to that of 11 [Eqs. (|l7|)-(|2^)] with the following two modifications. Since the irreducible polarizability does 
not include dynamic Coulomb screening (spin-density excitations are unscreened by the spin-conserving Coulomb 
interaction), we set J7i^,„„ = in Eq. (pl|). The second change concerns the vertex correction U^'^^^, which for 
spin-density excitations is given in the TDLDA by 

U^l^n = -/^^ / dz' M^)<Pj{^)^^{^)S{z ~ z') Pm{z')Mz') (23) 

instead of by Eq. (p^. In this equation, m{z) = n|(z) — f^x(•^) local spin density, n| and rii being the spin-up 

and spin-down local densities, respectively. 
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The CDE and SDE energies are given by the poles of fl and 11, respectively, which occur when the determinant 
l^ijUij^mn — Sij,mn\ vanishcs. In the numerical work presented in this paper, we solve this equation keeping all the 
subband levels. If, as in the quantum well structures we consider, the lowest two subbands are well separated in 
energy from the higher subbands, then one can approximate this determinental equation in the limit of low densities 
and temperatures by keeping only subbands 1 and 2, yielding 

(n02 + n°i)c/i2.i2 = 1. (24) 

For q ^ 0, this condition gives the resonance energies in the familiar form of Ando:0 

h^Q^i = e\ + 2e2ii7i2,i2(Afi - iVz), (25) 
where £21 = £2 ^ £1 and C/12,12 = ^^12.12(0) ~ ^viii the case of the CDE, and C/12,12 = ^V^2\2 foi" the SDEs. 



2. Polarized Electron Gas (LSD A ) 



In this Subsection, we introduce a generalization of the local-density approximation to density-functional theory 
discussed in Sec. [lAl which allows for different populations of the two spin orientations, i.e., a finite spin polar- 



ization. This local-spin-density approximation (LSDA) formalism is also based on the self-consistent solution of the 
Schrodinger-like Kohn-Sham equation, coupled with the Poisson equation and a local exchange-correlation potential. 
The main technical difference between LSDA and LDA is that the effective exchange-correlation potential in LSDA 
depends on the local spin polarization as well as the electron density. Therefore, one has to solve two Kohn-Sham 
equations, which contain spin-dependent effective potentials—for the two components oLthe spinor wave function. The 
LSDA was first formally justified by von Barth and HedinEj and Pant aiid RajagopaM and is suitable for studying 
ferromagnetic systems either with or without an external magnetic field.Eil 

To put the LSDA approach in context, we briefly review the theoretical and numerical evidence that the uniform 
electron gas in two and three dimensions embedded in a uniform positive background (the jellium model) undergoes 
a ferromagnetic transition at a certain critical density. A simple theoretical estimate for the density at which a 
ferromagnetic state will form may be obtained from Hartree-Fock theory, which treats the electron-electron Coulomb 
interaction to first order. For a uniform electron gas in three dimensions with spin-up and N- spin-down electrons. 



the ground-state energy in this approximation-ean be written in terms of the total number of particles N 
and the magnetization ra = — N^)/N aa^j 



N+ + N. 



Ne^ 3 /97r 
"2^10 \ 2 

Ne^ 3 /97r 



2/3 



1 



2/3 



1 



5/3 



4/3 



1 



5/3' 



4/3- 



(26) 



where ao = h'^/m*e'^ is the Bohr radius and = {SV/AirNy^^ /ao parameterizes the density. The first term in 
this expression is the kinetic energy, which prefers the paramagnetic state, while the second term is the exchange 
energy, which prefers to polarize the spins. At densities satisfying > 5.45, the exchange energy dominates, and the 
ferromagnetic state is stable. In two dimensions, a similar analysis yieldaiS 



p2_D 



2an 



1 



4V2 



(1 



^3/2 



V3/2 



(27) 



where now EE {A/TTN)^^^/ao. In this case, the condition for a polarized ground-state, Ejj^{rs,m = 1) < 
Ejjp{rs,m — 0), is satisfied if > 2.01. Thus, within the Hartree-Fock approximation, the spin-polarized state 
occurs at higher density (lower rg) than in three dimensions. 

Hartree-Fock theory neglects contributions to the energy beyond the exchange term and is therefore expected to 
overestimate the density at which the ferromagnetic transition occurs. Including these so-called correlation terms 
can only be done approximately, however. Currently, the most accurate method for performing these calculations 
are numerically intensive Monte Carlo techniques. CeperleycJ calculated the ground-state energy of an electron 
gas in two and three dimensions employing variational Monte Carlo (VMC). He found that in both two and three 
dimensions there is an intermediate density regime where a fully polarized state has the lowest energy compared to 
the unpolarized quantum liquid and the Wigner crystal. In 3D, the polarized phase is stable for 26 < < 67, and, in 
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2D, for 13 < Ts < 33. Additional results seem to indicate that in 3D there is a transition to a partially polarized liquid 
at Ts ~ 20 and to a fully polarized phase at ~ 50. c3 Eleven years later, Tanatar and CeperleyCj recalculated the 
ground-state properties of the electron gas in two dimensions employing the VMC technique and the more accurate 
fixed-node Green's-function Monte Carlo (GFMC) technique. The VMC technique predicted again a transition from 
the unpolarized to the polarized liquid at Ts between 10 and 20, consistent with Ceperley's results. The more accurate 
GFMC technique predicted a transition from the unpolarized liquid to the Wigner crystal at — 37, without an 
intermediate polarized phase. However, the authors point out that, near the transition, the polarized phase has an 
energy very close to the energy of the other phases and that, due to finite size effects and errors associated with their 
approximation method, their conclusion should not be taken as definite. This leaves open the possibility of a stable, 
fully-polarized phase in the two-dimensional electron gas. 

The LSDA falls somewhere between elementary Hartree-Fock theory and Monte Carlo calculations in terms of the 
quantitative accuracy of its predictions. Its main strength, and the reason we use this technique here, is that it 
can describe the inhomogeneous electron gas which exists inside quantum well structures and can therefore provide 
estimates of the critical density of the ferromagnetic transitions in these structures that could guide future experiments. 
In addition, the LSDA is a direct extension of the LDA and so allows a comparison between the two calculations 
which will be important in ruling out ferromagnetism as the source of the spin-density-excitation softening which 



appears in our TDLDA calculations [Cf. Sec. III]. We note that a similar problem to that of the spin polarization of 
the ground-state in quantum wells is the problem of "valley condensation" in Si-Si02 systems, where, instead of spin 
states, the electrons can occupy different valleys oLthe Brillouin zone. This problem has been studied in the past 
with techniques similar to the ones employed here.E3 In addition, the LSDA has been employed to study spin effects 
in wide parabolic quantum wells in the presence of a perpendicular magnetic field.cEl 

Our computations in the LSDA proceed as follows. After factorizing the complete single-electron wave function as 
was done in Eq. (^, we write down the z-dependent Kohn-Sham equation: 

'^'^ +Vc{.z) + VH{z) + V^c{^)]<t>l{z)^eMz), (28) 



2m* dz'^ 



where n is the subband index and a denotes the spin orientation, which can be up or down (-). We assume that 
no more than two subbands are populated and let n = 1 denote the lower-energy symmetric level (S) and n = 2 the 
higher-energy antisymmetric level (AS). We therefore need to consider four wave functions {4>i , 4>i, ^J, 02^) and 
their corresponding energies in Eq. ( p^ ) . 

This equation also contains the exchange-correlation potential, which in the LSDA formalism depends on both the 
density n{z) and the spin polarization m(z), which is defined as 

^^,^^nti-) + nti^)-n^i^)-n^i^)^ (29) 
n{z) 

In our calculations, we use the jjLarametrization of the exchange-correlation energy for the uniform 3D electron gas 
obtained by Ceperley and Alder:Ej 

where i = U (unpolarized, m = 0) or i = P (polarized, m — 1). The exchange-correlation contribution to the chemical 
potential is 



rs ^ (l + /31V^-f /3^r,)2 ^ ' 

The parameters in the previous expressions, as obtained by Ceperley and Alder,0 are: — —0.9163, = —1.1540, 
d'^ = -1.2218, d^ -1.5393, 7^ -0.1423, 7^ = -0.0843, /3f = 1.0529, P[ = 1.3981, = 0.3334, and 
= 0.261-L. For intermediate polarizations (0 < m < 1), we use an interpolation formula proposed by von Barth 
and IIedin,L3 in which the correlation energy has the same polarization dependence as the exchange energy: 

excirs,m) = exc(^^) + /(^") i^xcirs) - exc(^^)) (32) 

and 
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where 



V^c{r,,m) = V^cirs) + f{m) {V^cir^) ~ V'ic^r,)) + 

(4c(''.) - 4c(^s)) (sgn(^) - ^, (33) 



The Hartree potential Vh(z) is calculated as in the unpolarizcd formalism described in Sec. II A; it satisfies Poisson 
equation, Eq. (||), and is given by Eq. (^. 

To complete the specification of the problem, we note that the density associated with each subband and spin 
orientation is given by 

(35) 

where is the occupancy of each level, which at zero temperature is given by 

nl = N^{eF-el)Q{eF-el). (36) 
The total electron density may be written 

n(z)=^<(z), (37) 

n £T 

which implies that the Fermi level Ep is implicitly determined by the condition 



/oc 



[cf. Eqs. (pl)-(p^)]. The self-consistent solution of these equations proceeds exactly as in the LDA case with the 
difference that now, in each iteration, one has to solve two Kohn-Sham equations, Eq. (|28|), for the two spinor 
components of the wave functions, 05^. 

B. Self-Consistent Hartree-Fock Theory 

While conventional density-functional theory is a fairly accurate method for determining the properties of semi- 
conductor hetcrostructurcs in their paramagnetic and ferromagnetic phases, it is unable to address more complicated 
magnetic ordering such as antiferromagnetism. To study such phases, it is useful to return to a model Hamiltonian 
of the electronic system and search for the existence of broken-symmetry states within a mean-fiftld theory. In some 
cases, such as superconductivity, the mean-field theory gives a quantitative account of these phases.cZI More commonly, 
however, it sacrifices quantitative accuracy in favor of qualitative insight. This insight manifests itself not only in a 
physical intuition about the nature of the new ground state but also in the ability to study the distinctive features 
of the broken-symmetry phase, which may serve as a guide for interpreting experimental data in these systems. We 
adopt this point of view in what follows. 

To that end, consider a three-dimensional electron gas interacting through a potential V^(R — R') and confined along 
the z-direction by a potential Vconf{z)- A confining potential of this type is shown in Fig. |l| for a double quantum 
well, but the precise shape is unimportant for the development of the formalism. Given a particular Yconf{,z)^ one 
can construct its eigenfunctions i,n{z) and eigenenergies e„, n = 1,2, 3..., by solving the time-independent Schrodinger 
equation 



+ Vconf{z) 



^niz) = en^niz). (39) 



2m* dz^ 

In terms of these eigenfunctions, the quasiparticlc annihilation operator ip^CR) can be written as 



^.(R) = ^ E e''"""e.(^)c«k., (40) 
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where R = (r, z) = (x, y, z), k = {k^, ky), A is the transverse area of the sample, and c„kCT annihilates a quasiparticle 
in subband n, of transverse wave vector k, and with spin projection a (these conventions will be used throughout this 
paper). Note that, unlike the density-functional approach, the total confining potential Vconf{z) is specified at the 
outset and is not determined self-consistently. 

Defining a composite subband and spin index a — (rtg, gg) with summation over repeated indices implied, the 
Hamiltonian may be written in the basis defined by Eq. as 

H = Ho + Hint 

= ^ f^ak clk^ak + ^ ^ Kd,bc(q) clk+qCfek'-qCck'Cdk- (41) 
k kk'q 



In this expression, the quasiparticle energy 



Cak = e« + ~ M (42) 



is measured with respect to the chemical potential /i, and the matrix elements of the interaction are 

__ ^.q.(r-r') (^)^^^ ^(j^ _ j^,) (^,)^^^ (^,) ^ (43) 

Our goal is to solve this model within mean-field theory allowing for the possibility of broken-symmetry phases. For 
reasons that will be discussed in Sec. 0, we shall restrict attention to those ground states which are translationally 
invariant transverse to the layering direction, but we shall allow for off-diagonal order in both the subband and spin 
indices. This assumption excludes from the outset the study of intra-well charge- or spin-density waves and Wigner 
crystallization, and it is also implicit in the density-functional calculations discussed in the preceding Subsection. One 
could modify our treatment to include such phases, but the present model is sufficient for the purposes of exploring 
the effects of the inter- well degrees of freedom. The assumption of translational invariance implies conservation of the 
transverse wave vector, and so the quasiparticle propagator can be written 

/J 



Gabikn) = - dr e-"- (r. [c,k(T)4k(0)] ) , (44) 



where fc„ = (k, io-in), (3 — 1/T (h — ks — 1 throughout this paper), and the rest of the notation is standardcS 

The mean-field theory for our model is constructed by using this propagator to compute the electronic self-energy 
in the self-consistent Hartree-Fock approximation. This approximation corresponds to expanding the self-energy to 
one-loop order in the interaction and is shown diagrammatically in Fig. ||(a). The resulting self-energy is 

S„fc(k) = jYl KtMO) - VacMbO^ ~ k')] G,d{k'^). (45) 

The self-consistency of this approximation arises because the propagators used in Eq. (^5|) are dressed by the same 
self energy according to the Dyson equation 

[{iUn - Cak) Sab ^ ^ab{^)] Gbcihi) = Sac- (46) 

To completely specify the system of equations, the chemical potential is determined from the band-filling constraint 

Ns^jY. e-'^'-°- Gaaikm). (47) 

An alternative form of these equations, which will turn out to be convenient for future work, is obtained by inverting 
the Dyson equation for the interacting propagator Gab{kn) [Eq. (^)] by making an appropriate choice of the basis. 
From Eqs. ( ^3| ) and (^5|), the matrix Ca^Sab + 5^ab(k) is Hermitian and so possesses a complete and orthonormal set 
of eigenfunctions (pjj(k). These eigenfunctions satisfy the eigenvalue equation 

[eauSab + ^ab (k)] ft (k) = (k) (^^ (k) (48) 

(no sum on c) and the orthonormality relations 
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</.^(k)(^r(k) = Sab (49) 

and 

^^(kv:''(k) = 5'^^ (50) 

Eq. (H) is diagonal in the basis of these eigenfunctions by construction, so we may invert the equation to obtain 

lUJn - £'''(k) 

Substituting this relation back into the equation for the self-energy [Eq. (^5|)] and performing the sum over Matsubara 
frequencies yields 

i],,(k) ^jYl [y-bAcio) - VacM^ - k')] ^^(k')¥>r (k') /(i?'=(k')), (52) 

k' 

where f{x) = l/{e^^ + 1) is the Fermi function. The band filling constraint Eq. ( p7[ ) may similarly be written 

^^ = lE/(^'(k))^^n^ (53) 

ck c 



Eqs. (}48|)-(|53|) are the equations we will ultimately solve for a simple model interaction in Sec. VA. 

The eigenfunctions <Pa(k) and eigenenergies i?^(k) of the operator eakSab + ^^^(k) are the wave functions and 
energies of the quasiparticles of the interacting system within self-consistent Hartree-Fock theory. An alternative way 
of saying the same thing is that we have performed a mean-field decomposition of the Hamiltonian and diagonalized 
the result with a Bogoliubov transformation in the particle-particle channel. The annihilation operators for the 
interacting quasiparticles 7^ are therefore obtained from the bare operators via 7^ — ip'^(k)cak- 

In addition to the physical insight afforded by rewriting the self-energy equations in terms of these functions, 
the calculation of the energy, entropy, and specific heat of the system becomes straightforward. The energy is the 
expectation value of the Hamiltonian [Eq. ([4l|)] with the energy shift due to the chemical potential removed: 

E={H + fiN) 

^ {Ho) + {Hint) + fiN,. (54) 

Within our self-consistent Hartree-Fock theory, the contribution to the energy from the interaction term in the 
Hamiltonian is shown graphically by the diagrams in Fig.|(b). These diagrams lead to the result 

(^int) - ^ E ^baikn) ^abC^), (55) 



2 



which allows us to write Eq. (^^ as 



E^^J2 e"*""°" [2e„k'5ab + S,b(k)] GbaiK) + M^s- (56) 



2 

k„ 



Substituting Eq. (|5l|) into this equation, using the eigenequation Eq. ( ^8| ) and orthonormality relation Eq. (pO|), and 
performing the Matsubara sum, the total energy becomes 

E = IY1 [l'^a(k)|' eak + E^ik)] fiE^ik)) + f,N,. (57) 

k 

The entropy in the interacting basis is simply the standard free-fermion result 

S = -J2 ifiE^k)) In fiE^k)) + [1 - fiE%k))] In [1 - /{E^k))]} , (58) 

k 

and the specific heat is obtained directly from this equation: 
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dT 

t 

In order to get a full picture of the interacting system, one must go beyond the single-particle properti es and 
thermodynamic functions and examine the response of the system to external perturbation. As discussed in Sec. II A 1 , 
resonant inelastic light scattering has proven to be a powerful tool for studying the charge- and spin-density excitations 
in semiconductor heterostructures. Consequently, we will focus on the generalized density response function and the 
resulting collective excitations which can be observed in these experiments. As in Sec. [I A 1, we will not compute the 
form factors necessary to connect the polarizability to the inelastic light scattering cross section; we merely note that 
the collective modes we will discuss are detectable in these experiments with a particular (and structure-dependent) 
arrangement of scattering angles, polarizations, etc. 

We begin by defining a generalized density operator 

p'^(R) = 7^t^(R)<,,^.,(R), (60) 

where = {cr°, cr^, cr^, cr^} ~ {1, cr^, cr^, cr^}. Note that the number density ri(R) = p''(R) and the spin density 
s*(R) = (?i/2) /o*(R), i = 1, 2, or 3. Suppose we can couple to this density through an external force i^gxt(R-7 which 
adds a term 

iJext(i) = dRp^(R,i)^cxt(R,i) (61) 

to the Hamiltapian (a sum on /i is implied). The linear response of the generalized density to this perturbation is 
then given byCJ 

(5p''(R,t)) = % J dR'dt'eit') ([p''(R,i),p'^(R',t')]) i^cxt(R',i'), (62) 

where the angle brackets denotes the thermodynamic average in the absence of -ffoxt (t) . 

As mentioned above, we assume that the interacting system is translationally invariant in the transverse direction, 
and it is also time-translation invariant. Thus, we may introduce partial Fourier transforms 



p'^(R,i) = l^ / ^e'(''--"*V^(r,z,t) 



(63) 



and write Eq. ( |62| ) as 

{Sp'^{<l,z,u;))^i I dz'dt'e'^'' e{t') {[p^{ci,z,t), p''(-q, z', i')]) F^^.ici, z' , uj) , (64) 



From Eqs. ( [40| ) and (pOD, we can rewrite the Fourier components of the density operator as 

P^{^,z)=CnAn,<a,PaM) (65) 

with 

Pafc(q) = Xl"^akC6k-Hq- (66) 
k 

Using this relation, we obtain the final form for the generalized density response: 

(5p^,(q, u)) = -n,fc,,,(q, u:) (q, u:) (67) 

where 



na6xd(q, ^) = y dt e'"* ( [p,,(q, t), pL(q. 0)] ) , (68) 
(<5p^(q,z,^)) =C(^)C..(^X,.. {5paM.^)), (69) 
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and 



F,r(q,^) = J dz'^,Jz')CAz')<^a, ^^ext(q,^',^)- (70) 

In order to obtain the subband- and spin- resolved polarizability [Eq. (|6^)], we compute the polarizabihty in Mat- 
subara frequencies zf„ = 27rnT, 

HabMln) = - t dre^""" (Tr kfc(q, r)pt^(q, 0)1 ) (71) 







{in = (q, ii^n)), and analyticaUy continue to real frequencies by the convjeational substitution ivn — » + i5 (ci. 
Ref. kg). This polarizability is calculated within a conserving approximationEj using the diagrams shown in Figs. y(c) 
and 0(d). These diagrams yield the expressions 

+ Qn) (72) 



A- 

k„ 



for the polarizability and 



lab,cd{krn, km + (J„) = SacSbd ^ ^ l^bf.eaO^ - k') - Vf,a,e/(q)] 

K 

X Ggeik'i) Gfh{k'i + qnhghMk'l, k'l + qn) (73) 

for the vertex function. Contrary to the usual convention, we include the RPA screening diagrams in the vertex 
function and thus do not distinguish between reducible and irreducible polarizabilities. The reason is that, in the 
broken-symmetry phases we will examine, the propagators are no longer diagonal in spin, so these screening terms 
enter into both the spin- and charge-density response functions, reducing the usefulness of the division between 
reducible and irreducible polarizabilities. Note that even these approximate equations for the polarizability must be 
solved numerically for a general interaction. If the interaction is wave- vector-independent, on the other hand, the 
vertex equation can be solved simply (cf. Sec. ^). 

The interpretation of the polarizability as a response function can be used to visualize the real-space density 
fluctuations it describes. If one perturbs the system with the external force Fpxt(R-i^) °^ e*^'^''"'^'^ Eqs. (|6^)-([70|) 
imply that the response in real space will be 

(<5p^(R,i)) « e'^^'i-'— *)C„(^)e„,(^X., n,,^,,(q,c^)/°f , (74) 

where /^^* is related to how this force perturbs the subband and spin indices. For example, to examine the spin-density 
response between the lowest two subbands, 

flf = ['5nei5„,2 + '5„^2'5„,i] 5„^i^5„^i. (75) 

This procedure works for general wave vector and frequency, but has special significance when these quantities 
correspond to a collective mode of the system, in which case the response gives the polarization of the collective mode. 
Note that, since the response function diverges at a collective excitation, in practice one applies this formula by adding 
a small imaginary part to the frequency to control this divergence. As with the equations for the polarizability, this 
numerically intensive approach for obtaining the collective mode polarization is required for the broken-symmetry 
phases where the off-diagonal terms make an analytic calculation difhcult. 

III. INTERSUBBAND INSTABILITIES 

In this Section, we shall study the instabilities of an electron gas confined to coupled double-quantum-well (DQW) 
and wide single-quantum-well s truc tures by examining their intersubband collective excitations. We begin by applying 
the formalism described in Sec. [I A to the typical GaAs / AlxGai_xAs DQW structure shown in Fig. |l|; similar samples 



were used in the experimental study of Ref. Although our TDLDA calculation includes all the subband energies 
shown in the figure, we are only concerned with intersubband transitions between the lowest two subbands (subbands 
1 and 2), whose typical separation /S.sas = £2 — £i ~ 1 meV. The exact values of IS.sas are plotted in Fig. ||, along 
with the calculated long- wavelength energies of the intersubband collective modes as a function of the sheet density. 
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Two distinctive features of these collective modes are seen in Fig. I The first feature is that the CDE energy 
becomes less than Ag^s around Ns « 0.2 x 10^^ cm~^rx-| This behavior, originally predicted in Ref. |l^ for single 
quantum wells, has recently been verified experimentallyE3 and will not be discussed further. The second - and for 
our purposes more important ~ feature is that the intersubband SDE energy goes to zero below a critical density 
Nc ~ 0.7 X 10^^ cm^^ and becomes finite again below 0.1 x 10^^ cm^^. This softening of the intersubband SDE 
mode indicates that the normal system with a "metallic" Fermi surface is unstable at or below Nc, since it can 
spontaneously create spin-reversed intersubband electron- hole pairs ("triplet excitons") at no cost in energy. We 
conclude that there is a phase transition in this DQW at the critical sheet density Nc from the normal 2D electron 
liquid to a triplet intersubband exciton liquid and that the system re-enters the normal phase at a lower density. 
This electronic phase transition is due exclusively to exchange-correlation effects which make the vertex correction or 
excitonic shift larger than As^s [Eq. (p5[)] and thus cause the SDE to collapse. The re-entrance of the normal phase 
at lower density is simply explained by the fact that the vertex correction vanishes as iV^ ^ according to Eq. (|2^). 
We note that the transition to a Wigner crystjalraccurs at much lower Ng values than those considered in this paper 
and so does not account for the SDE collapse. liacJ 

Additional evidence in support of the excitonic transition comes from the density dependence of the mean-field vertex 
correction | UxcXi2 Ij which is plotted in Fig. 3(a). This vertex correction consists of the spin-polarized, exchange- 
correlation-induced vertex function Uxc [Fig. |(b)] and the uncorrected intersubband polarizability Xi2 = +n2i 
[Fig. ^(a)]. In the two-subband limit, the vertex-corrected irreducible response function is given by Xi2[l~t^xcXi2]~^j 
which clearly has an instability when | UxcXi2 1^ 1 [<^f- Eq. (|2j)]. As shown in Fig. ^(a), this "Stoner criterion" is 
satisfied in the range of sheet densities in which the SDE has collapsed. 

Having established the existence of the SDE-coUapsed phase for a particular DQW structure, we can vary the 
geometry of the structure to study the persistence of this phase. For well widths fixed at 139 A, the calculated 
zero-temperature phase diagram in terms of the sheet density and barrier width ds is shown in Fig. ^. For very 
small barrier widths, Ag^g is too large for the vertex correction to overcome it, even at low densities, making the 
normal phase the only stable phase. For very large barrier widths, on the other hand, Asas is exponentially small, 
and the normal phase gives way to the excitonic one at extremely low densities. Of course, for such exponentially 
small AsaSj the critical temperature for our predicted instability is also exponentially low, implying that the phase 
transition in large-ds, low-density DQW structures would be unobservable in practice.EJ At higher densities in the 
large-ds structures, the normal phase re-asserts itself due to the diminishing influence of the vertex correction. For 
intermediate values of the barrier width, we see the re-entrant behavior described above with an upper critical density 
that is nearly independent of the barrier width. A similar phase diagram is obtained when the barrier width is fixed 
at 40 Aand the well widths are varied, as demonstrated by the inset to Fig. By expressing the sheet density and 
AsAS in dimensionless form, the data from Fig. |^ can be reduced to the universal zero-temperature phase diagram 
shown in Fig. |^. From this figure, we see quite clearly that the excitonic phase appears in a range of densities below 
a critical symmetric-antisymmetric splitting Ag^s and that this density range increases as Ag^s decreases. 

The results discussed so far for the intersubband SDE energy are for zero in-plane momentum transfer. Fig. |^ 
shows the dispersion relation of the lowest intersubband spin-density mode for several sheet densities Ns above Nc- 
The most important feature shown is that the SDE mode becomes soft at a finite value of the in-plane momentum 
transfer, Qc, at a critical density somewhat higher than the critical density found at zero in-plane momentum transfer. 
The critical momentum is given hy Qc = kpi — kp2, where kpi is the in-plane Fermi wave vector of the ith subband in 
the normal state. This indicates that the excitonic instability, for densities slightly above Nc, may occur at a finite 



value of in-plane q, a point to which we shall return in Sees. VB and VI 



Thus, we expect that a softening of the intersubband SDE may be observed either at fixed q and varying Ng [cf. 
Fig. |^(a)] or for q ^ Qc at certain fixed values of Ns [cf. Fig. |^(b)]. In Fig. ^(a), we show the spectral function of the 
intersubband spin-density excitations as the instability is approached from the high-density side at very small q. As Ns 
approaches Nc, the spectral peak shifts towards zero energy and the line narrows. On the other hand, the approach 
to the critical momentum transfer qc at densities just above Nc, presented in Fig. ^(b) for Ns — Nc = 0.7 x 10^^ 
cm~^ at various values of q approaching qc, the line broadens as it softens. Thcs p B Pt j ctra, are similar to what would 
be observed in Raman scattering experiments in the cross-polarization geometry.cSEjOlEJ For sheet densities within 
the SDE-collapse d pha se, the TDLDA treatment employed in the Section breaks down, and we must employ the 
techniques of Sec. [IB to determine the collective modes, which we shall do in Sec. 

Finally, we mention that the instability discussed here for DQWs can also happen in wide single quantum wells, 
in which the effective potential develops a barrier in the center at moderate densities, becoming similar to a DQW. El 
An example of this phenomenon is illustrated in Fig. ^; notice that the electron density profile is similar to that for 
the DQW system shown in Fig. |l]. In this effective double well system, the separation between the two lowest lying 
subbands is small, about 0.2 meV, as in the coupled DQWs, and the possibility of a suppression of the intersubband 
SDEs by the vertex correction arises. A calculation of the SDE energies shows that this is indeed the case, as can be 
seen in Fig. For the square well considered, whose width is 1000 A, the region in which the SDE has collapsed 
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is somewhat smaller than that for the DQW of Fig. |l| with a 40 Abarrier width, but it occurs around the same 
Ns ~ 0.4 X 10"'^^ cm~^. Hence, inelastic light scattering experiments should be able to detect the SDE instability in 
wide square wells also. 



IV. FERROMAGNETIC PHASE 



In the preceding Section, we showed that the electron gas confined in a coupled double-quantum-well structure has 
an instability which is indicated by the vanishing of the energy of the intersubband spin-density excitations. One of 
the goals of this Section is to establish whether this instability corresponds to a phase transition from the normal 
spin-unpolarized (param agnet ic) ground state to a spin-polarized (ferromagnetic) one. A second goal is to employ the 
LSDA formalism of Sec. [I A 2 to the ferromagnetic phase transition in single square quantum wells. In particular, we 
are interested in the dependence of the critical density on the well width, as it would reveal an aspect of the crossover 
from two- to three-dimensional behavior of the electron gas which has not been explored in the past. Our main finding 
is that the LSDA does predict a ferromagnetic transition in single quantum wells and, moreover, that as the electron 
gas is widened the transition density decreases (the ferromagnetic phase becomes less favorable). This dependence 
of the critical density on the electron-layer width agrees with the well-known fact that the influence of the Coulomb 
interaction is stronger for lower dimensionality. 

We first determine whether the excitonic instability studied in Sec. [II can be explained in terms of a ferromagnetic 
transition, that is, a transition from the usual spin-unpolarized ground state to a partially or fully polarized one. We 
concentrate on a coupled double-quantum- well system in which each individual square well is 139 Awide and the well 
separation is 40 A~ the main structure studied in Sec. [II. As in the better-known 3D case, a fully polarized phase 



is expected at low density, and a normal, unpolarized phase at high density. To compute the spin polarization of the 
system at a given sheet density within the iterative, self-consistent LSDA method, it should be enough, in principle, 
to introduce a slight asymmetry in the initial choice of spin densities. If the correct ground state were unpolarized, 
the initial asymmetry would rapidly disappear in the iteration process. On the other hand, if a polarized state were 
expected, the initial small polarization would increase until convergence to the fully polarized is achieved. However, 
since the solution of the self-consistent set of equations of the LSDA method is affected by numerical inaccuracies, in 
practice our algorithm is sensitive to the initial guess for the spin density profiles. Therefore, we employ the following 
method to determine the ground-state polarization of the system. For a given density iV^, we solve the self-consistent 
algorithm starting with the spin densities ni{z) = rinu{z), ni{z) = (1 — ri)nu{z), nj(z) = 0, and n^iz) = for 
rj = 0.55 and 0.95 and n„(z) the density from the unpolarized LDA calculation (only one subband is occupied at the 
low Ns studied). If the calculation converges to a polarized or an unpolarized state for both values of ?7, we take the 
result as the true polarization at that given Ns- Otherwise, we assume that the result is affected by the insufficient 
numerical precision. 

With this method, we obtain the phase diagram of spin polarization as a function of Ng for the coupled DQW 
defined above. The calculation converges to a fully polarized state, regardless of the initial choice of spin densities, 
for iVg < 4 X 10® cm~^ and to an unpolarized state for Ng > 1.7 x 10^ cm~^. In the range Ns = 0.4-1.7 xlO^ cm~^, 
the solution is polarized for the choice rj = 0.95, and unpolarized for rj — 0.55. 

Therefore, our first conclusion is that, within the LSDA, there is a ferromagnetic transition as a function of the 
electron density Ns in a coupled DQW system; for our particular choice of parameters, the transition occurs around 
Ns ^ Ix 10^ cm^^ at zero temperature. On the other hand, the ferromagnetic transition occurs at a density which is 



almost two orders of magnitude smaller than the critical density Nc of the excitonic transition of Sec. III. Moreover, it 



does not seem to show re-entrant behavior at lower density as the excitonic transition does. Based on these differences. 



we conclude that the ferromagnetic transition cannot be identified as the excitonic phase transition of Sec. [II 



Having achieved the first aim of this Section, we now turn to the second, namely the possibility of ferromagnetism 
in single square wells. We wish to determine the critical density of the ferromagnetic transition as a function of the 
well width in order to take the first step in studying the evolution of this transition as a functio n of th e dimensionality 
of the electron gas. The analytical results for two and three dimensions summarized in Sec. 11 A 2 indicate that, in 
2D, the exchange energy is more important than in 3D, indicating a higher 2D critical density. In quasi-2D systems, 
we therefore expect that the ferromagnetic critical density should decrease as the well width is increased. 

Employing the method described above, we study the ferromagnetic transition in five single square wells of widths 
dw = flOi 2ao, 4ao, Gag, and Soq, where oq — 98.7 Ais the effective Bohr radius in GaAs. The resulting phase 
diagram is presented in Fig. ^ in terms of Ns versus well width, and also in terms of the 2D and 3D and well 
width; the vertical bars give the density range where the polarization of our solution depends on the initial choice 
of spin densities. As expected, we see that the critical density decreases with increasing well width. However, the 
limiting values of rs for narrow and wide wells cannot be directly compared with the pure two- and three-dimensional 
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Hartree-Fock values for three reasons. First, the analytical results are obtained in the jellium model, whereas in our 
quantum well calculations the positive charges of the ionized donors are located far away from the electron gas, which 
should produce an important change in the direct Coulomb energy and hence affect the ferromagnetic critical density. 
Second, our calculation includes correlation effects which go beyond the exchange interaction. And the third reason 
to expect differences is that, in the LSDA, exchange is treated in a local and static approximation. 

This exchange-correlation potential is obtained from the ground-state energy of a uniform electron gas, which is not 
known exactly and must be calculated within some approximation scheme. Several authors have proposed different 
parametrizations of the exchange-correlation potential, which are based on different calculations of the 3D-electron- 
gas ground-state energy and which therefore differ somewhat from each other. We have checked some of our results 
using the parametrization of Gunnarsson and Lundqvist,E3 which gives a larger difference between the potentials 
for different spin polarizations than the Ceperley- Alder one, making the ferromagnetic phase more favorable. For 
example, for the square well of width dw = 4ao, the critical density is 1.3 x 10^" cm~^, ten times higher than with 
Ceperley- Alder. Thus, our calculations indicate that within the LSDA a ferromagnetic transition is indeed present in 
the quasi-two-dimensional electron gas, but the value of the critical density for the transition cannot be accurately 
determined with this method. 



V. ANTIFERROMAGNETIC PHASE 



Based on the results of the preceding two subsections, we can draw the following conclusions: (1) double or wide 
single quantum well structures can undergo an excitonic instability in a range of densities, and (2) this instability 
cannot be associated with a ferromagnetic transition. We would like to understand the nature of the ground state 
of the excitonic phase and determine its excitations, but we see that the density functional approach cannot assist 
us further. In order to make progress, we employ the self-consistent Hartree-Fock formalism developed in Sec. II B, 
which allows us to study ground states with broken symmetries which are not accessible through ordinary density 
functional theory. 

As pointed out in Sec. II B , the self-consistent Hartree-Fock model for systems with a Coulomb interaction is not 
quantitatively accurate but should reproduce the qualitative features of the new phase. Our goal is therefore to 
construct a minimal model of the systems which exhibit the excitonic instability which is still tractable. We can make 
several observations about the excitonic instability which will guide us in developing this model. First, the instability 
in the LDA calculations occurs with a wave-vector independent interaction, as seen from the LDA expressions for the 
irreducible polarizability, Eqs. ( p^ ) and (p^). Second, the instability is signaled by a soft mode which is an excitation 
between the lowest two subbands and is accompanied by a spin flip. Thus, our minimal model will focus only on 
the lowest two subbands, allow for off-diagonal terms in the subband and spin indices, and use an interaction which 
in constant in wave-vector space. In the rest of this Section, we discuss the repilts which follow from this minimal 
model; a preliminary report of this work has already appeared in the literatures 



A. Point-Contact Model 

A q-indepent interaction has several simplifying consequences. First, as seen from Eq. (^Sf), the self-energy will 
also be wave- vector independent. From Eq. (|4^), this result implies that the eigenfunctions <Pa(k) are also wavevector 
independent and that the interacting eigenenergies may be written 

E^(k) = e^ + 0-,. (76) 

Because these eigenenergies are the energies of the interacting quasiparticles, the effect of the interaction is to 
shift and/or rearrange the non-interacting bands without destroying their parabolic dispersion. Consequently, the 
broken-symmetry rphases available within this model will not exhibit any energyr-sap, in contrast to what occurs in 
superconductivityEll and in other excitonic transitions which have been studied.!^ The reason for this difference is 
that the interacting quasiparticles are linear combinations of electrons and holes in the latter cases, whereas in this 
case they are linear combinations of electrons from different subbands in the former. 

Second, the simplest interaction which is q-independent is a delta-function in real space: V(Ii) = V6(TC). Inserting 
this form into Eq. (^) yields 

Vab,cd{(l) — Sa^at^a^aa yn^n,,,n^na (77) 

with 
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K„„„n.n, J dzCnSz)i^^{z)CnS^)^,Jz). (78) 

Making the standard assumption that the quantum weU structure is symmetric, the wave functions may be 

chosen to be real and parity eigenstates. From Eq. (|7^), these symmetries imply that the order of the subband indices 
in the matrix elements is irrelevant and that matrix elements of the form Vli,i2, ^12, 22) and cyclic permutations 
vanish. Thus, there are only three independent matrix elements: Vu = Vu^u, V22 — V22,22, and V12 = V\i.ti + 
permutations. Since we are interested in the minimal model exhibiting the excitonic instability, we will neglect V\\ 
and V22 in order to focus on the intersubband effects, leaving a model with two free parameters: the sheet density N ^ 
and the intersubband interaction V\i. 

Applying these approximations to Eq. (^) , we find that the self-energy reduces to 

This equation is supplemented by the band-filling constraint, Eq. (p3[), with the electron density in each interacting 
subband given by 



= lE/(^'W)- (80) 

k 



In addition, the integral equations determining the polarizability [Eqs. (|7^)-(|7^)] can be reduced to a matrix 
equation. The first step in this reduction is to observe that the vertex function 7 in Eq. (73) is a function only of q„ 
when the the interaction has no q-dependence and satisfies the equation 

lab,cd{qn) = 5ac5bd + [Vba.ef - Vbf^ea] ^^ef,ghi'l^) Igh^cdi^n) (81) 

where 

ni°id(9«) = 5 E Gca(fc™) GM{k^ + qn) (82) 

km. 

is the polarizability without vertex corrections. In the interacting system, H^^^^^ may be written through the use of 
Eq. (|!|) as 

nit.(9") = v*aviv>:I n^^(g„). (83) 



After analytically continuing to real frequencies, the function 'iY^ {q) is just Eq. (j2^) with = 1 and using the 
interacting energies [Eq. (u6|)]. Inserting this result into Eq. (ra) and rearranging yields 



Inverting this matrix equation gives the polarizability. When this equation cannot be inverted for a particular (q, cj), 
i.e. when the determinant of the matrix in brackets vanishes, the system exhibits a collective mode. The polarization 



of these collective modes can be determined by, for example, the method discussed in Sec. II B, 

Before proceeding to the discussion of the normal-state properties of our point-contact model, it is appropriate to 
discuss briefiy the question of how well this model approximates reality. The Hartree-Fock approximation is known 
to be a poor one in the interacting electron gas because it neglects screening effects, but including these effects 
realistically is a difficult problem which has not yet been resolved. In computing the collective modes, this problem is 
amplified by the distinction one should draw between the interaction between the bubbles, which is unscreened, and 
the interaction within the bubbles that gives the non-trivial part of the vertex equation, which should be screened. 
Summing a particular set of screening diagrams may reduce the error introduced into the self-energy, but would render 
the collective mode calculation completely intractable. Thus, we adopt a strong approximation to the actual Coulomb 
interaction, a point-contact interaction, which should reproduce the qualitative features of the screening effects while 
leaving a solvable set of equations. 



16 



B. Normal-State Instability 



Our goal is to use the point-contact interaction model we have just described as the minimal model for the excitonic 
instability. It is therefore critical to verify that this model exhibits the instability in its normal state, an exercise which 
will also provide guidance to the nature of the broken-symmetry state. The first step in this process is to compute the 
self-energy and chemical potential in the interacting system. In the normal paramagnetic phase of this model, we know 
that the wave functions of the subbands are not mixed and that the band structure consists of two parabolic subbands, 
each with degenerate spin-up and spin-down components. Consequently, the self-energy is diagonal in subband and 
spin indices, and the eigenvectors of Eq. (^ ) are given by (p^ = (5„ with the corresponding eigenvalues satisfying 
giT — gii ^ a,nd e^^ — e^^ = 62- It then follows that n^^ — n^^ = ni, n^^ — -n?^ = n2, = = ^i' 

and S2T,2T — ^2j,,2j, = ^2- Inserting these formulae into Eq. (|79| ) yields the self-energy equations Ei = Vi2n2 and 
S2 = Vi2ni. 

At finite temperatures, the simultaneous solution of these self-energy equations and the band-filling constraint 
Eq. (^) must be performed numerically, but at T = the solution can be obtained analytically. In order to exhibit 
the zero-temperature solutions in a compact form, we define the non-interacting subband splitting in terms of the 
eigenvalues of Eq. (|39|): Ag^g = £2 — ei. The interaction renormalizes this splitting without modifying the identity of 
the underlying subbands. Applying the self-energy equations allows us to write this renormalized splitting as 

AsAS = 62-61 

= A°s^s + Vi2{ni~n2). (85) 

If the chemical potential is measured from the bottom of the lower subband and fi < Ag^g, then only the lower 
subband is occupied at zero temperature, and we shall refer to this state as the A^i phase. Alternatively, if ^ > Ag^g, 
both subbands are occupied, and we shall refer to this as the A^2 phase. By solving the equations for Ag^g and the 



band filling constraint Eq. (53) simultaneously, we obtain 



/ A»»s + ijV,l/„ (JVi phase) 



and 




(Ni phase) 

f^^< "n" 1 a° (87) 
^ 1 Af. ■ 1 (iV2 phase) ^ ^ 

The crossover from A^i to N2 occurs when = Ag^g or Ns/2NoA'^j^g = 1/(1 ^ ^^0^12). 

With the self-energy and chemical potential computed, we are now in a position to study the collective intersubband 
spin-density excitations of this model to see if the excitonic instability occurs. The condition for these excitations 
obtained from Eq. (|4|) through the use of Eq. ( |8^ ) is 



1 ± Re [ni2 {q,u;)+U^\q,u;)]Vi2=0, 

where the upper (lower) sign is for the spin (charge) density excitation and we have dropped the spin indices since 
n°'' does not depend on them in the normal state. At q = Q-and T = 0, an analytic solution for the frequency of the 
collective modes cuq is possible which is just the Ando resultE3 with a renormalized subband splitting: 

^0 = (A^As)' T 2Fi2(ni - n2)A*As. (89) 

It is clear from this expression that the spin-density excitation (SDE) will soften completely when Ag^g < 2Vi2(ni — 
712). At zero temperature, the relations derived in the preceding paragraphs can be used to show that the softening 
occurs in the A^i phase when A^sVi2/2Ag^g > 1, and in the A''2 phase when A'oVi2 > ^. These boundaries are shown 
in Fig. |lj along with the line separating the A'^i and ^''2 phases. The shaded area represents the region in which the 
spin-density excitation is soft at g = in this model. The phase corresponding to the soft mode occupies a large 
region of parameter space and obtains at any density providing the interaction V12 is sufficiently strong. 

At non-zero wave vector, the E g. (|8^ ) must be solved numerically. From the low-density, A'^i phase, the results of 
this calculation are shown in Fig. n3^ for increasing density at fixed interaction V12. As illustrated by the figure, the 
q = spin-density excitation softens with increasing density and vanishes beyond a critical density which depends on 
the interaction strength V12. At higher densities, the small-g solutions to Eq. (p8|) become imaginary, and the largest 
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imaginary frequency - indicating the most unstable wave vector - continues to be at g = 0. Thus, the transition from 
the low-density side is very clearly a q = instability, as is also seen in the TDLDA calculations. 

The situation from the high-density, A'2 side is more complicated. The dispersion curves for this case computed 
from Eq. (^) at fixed density and increasing interaction strength closely resemble those in Fig. |7[ As the interaction 
strength V12 increases, the energy of the SDEs are reduced until the entire dispersion curve collapses along a range of 
wave vectors from q = to kpi — kp2, the difference of the Fermi wave vectors of the two subbands, at a critical V12. 
The reason for this collapse is seen from Eq. (pSh when one observes that at T = in the 7V2 phase 



Re [n^^ (^q^Lj)+ tfi (g, w)] = -2iVo (90) 

for q < kFi+kp2, indicating that if the energy of a SDE with one of these wave vectors vanishes, they all must vanish. 
At interaction strengths slightly larger than that required for a collapse, the solutions to Eq. ( ^ ) are imaginary in 
this range of wave vector with the largest imaginary frequency occurring at q = 0. These calculations suggest a 



q = transition from the high-density side, but the density-functional computations in Sec. [II find a softening at the 
non-zero wave vector q = kpi — kp2- The discrepancy arises from the fact that the calculations in Sec. Ill include all 
the subband levels, whereas the model in this Section contains only two. The inclusion of higher subbands, even if 
they are unoccupied, disrupt the cancellations leading to Eq. ( po|) and yield a g-dependent function which will select 



some wave vector. The consequences of a finite ordering vector will be discussed in Sec. VI but are beyond the scope 
of our simple model. We nonetheless expect that many of the qualitative insights from our model will apply to more 
comprehensive theories. 

One of the most important insights that can be gained by the study of the normal-state instability is a physical 
intuition about the nature of the ground state in the region of the the phase diagram in which the SDE has softened. 
In structural phase transitions, a particular phonon softens, and the polarization of that phonon determines the new 
structure. In the same way, the polarization of the soft spin-density excitation should reveal a great deal about the 
new ground state. This polarization can be computed by means of Eqs. (|7^)-([75|), which leads to a density response 
in the normal state given by 

(p^(R,t))(xe'('l-"— *) 6(^)6 W [nl2(g,^)-Hn2l(g,^)] (f-iy), (91) 

which is shown graphically in Fig. |lj(a). Notice that the response occurs in the spin channel and has the character of 
a traveling wave transverse to the layering direction. Also observe that, since the wave function of the lowest subband 
^i(z) is even in z due to the assumed symmetry of the quantum well structure, and the wave function of the next 
subband ^2(2) is odd, the overall response odd in z. Thus, the intersubband spin density excitation corresponds to a 
spin-density wave in which the spin density is antiferromagnetically correlated between the quantum wells. 

With this interpretation, the nature of the excitonic instability becomes obvious: the softening of the q = SDE 
corresponds to the formation of antiferromagnetic order between the quantum wells with no modulation of the spin 
density within a well. A profile of the resulting spin density is presented in the inset of Fig. Note that there are 
three degenerate SDEs in the normal phase which may be associated with spin polarizations along the three Cartesian 
directions. Although all three soften simultaneously, we expect that a particular spin direction will be selected by the 
system as in the ferromagnetic transition, but there is a priori no restriction of the direction of the spin polarization 
selected. 

To summarize the results of the application of our simple model to the normal state, we find a softening of the 
intersubband spin-density excitations at g = for a wide range of model parameters. These collective excitations 
produce a spin-density modulation which is antiferromagnetically correlated between the wells. Thus, the softening 
of the g = spin density modes indicate the formation of an antiferromagnetic phase in which the spins in both 
wells are polarized but in opposite directions. In the more realistic density-functional calculations of Sec. |lll, a q 7^ 
instability was indicated which would imply an additional intra- well modulation of the spin density (see Sec. VI). For 
simplicity, we focus on the q — Q phase. 



C. Antiferromagnetic Ground State and Thermodynamics 

With an understanding of the nature of the excitonic instability as an antiferromagnetic ordering, we are able to 
examine the new ground state and its properties within our model. We first note that both spin rotation invariance and 
parity are broken in the antiferromagnetic state, so we are studying a genuine symmetry-breaking phase transition. 
The first task encountered in studying such a transition is to identify the order parameter. For concreteness, consider 
the antiferromagnetic phase with the spin density aligned along the x direction. From Eqs. (|6^) and (|65[), the only 
non-zero component of the spin density is then 
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(92) 



Because the antiferromagnetism requires that (s^(z)) be odd in z, only the terms which are ofF-diagonal in subband 
index are non-zero in this phase. Defining 



rnin2 



E 



(93) 



and using Eq. (66), this imphcs that 



'2kT''lki 



(94) 



^Uiz)Uiz) [P12+P21] 



(95) 



P12 + P21 finite in the antiferromagnetic phase and vanishes in the normal phase, so it is the order parameter of the 
phase transition. From the spin-rotational invariance of the system, the direction in which the spin density can be 
polarized in the antiferromagnetic phase is arbitrary, so the general form of the order parameter is 



(96) 



With this identification, we can construct a form for the self-energy matrix which permits the off-diagonal expec- 
tation values we require yet preserves the other symmetries of the system. This ansatz may be writtencil 



Si Sod 

Si Sod 

Sod S2 

Sod S2 



(97) 



Inserting this form into the eigenvalue equation Eq. (48) yields the eigenvectors 



Va 



u -V 

w 

w u 

w u 



(columns correspond to the superscript) and eigenenergies e'^ = {e+,e+,e ,e } [cf. Eq. (|76|)] , where 



It 



^SAS + S2 ~ Si 

2D 



-,1/2 



(99) 



AO 

^SAS ■ 



S2 -hS 



(100) 



and 



D 



^SAS 



Si 



1/2 



(101) 



The eigenenergy spectrum indicates that the single-particle band structure in the antiferromagnetic phase consists of 

two sets of two degenerate parabolic subbands separated in energy by A_| = e"*" — e~ = 2D. This band structure 

is similar to that in the normal phases, but the wave functions of the interacting quasiparticles are considerably 
different: from Eq. (^8|), we see that these quasiparticles are linear combinations of electrons from the two non- 
interacting subbands. 

The parameters in the eigenenergies and eigenvalues are determined from the self-consistency equations obtained 
by substituting the eigenvectors into reduced self-energy equation, Eq. (|7^). This procedure gives 
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Si = Vi2 (n-u^ + n+v^), 
E2 = V12 (n^u^ + rt^ti^) , and 
Sod — 2Vi2 {n^ — n+) uv. 

Using the definitions of u and and v [Eq. (p9[)] and the band filling constraint 

iV^ = 2 (n" + n+) , 

these self-energy equations can be written 



(102) 
(103) 
(104) 

(105) 



^SAS + S2 - Si - 



^SAS 



l~Vi2in- -n+)/2D 
Vi2{n~ — n+) 



and 



D 



od- 



(106) 
(107) 
(108) 



In the antiferromagnetic phase, Sod 7^ 0, implying from Eq. ( |108D that D — V1 2 (n — n+). Eq. 
^SAS + S2 - Si = 2AgAs. from which we obtain = (AgAsr + ^^od by Eq. ( |lOl| ). Comparing 
for D yields 



(107) then becomes 
the two expressions 



Sod - [V^il(n--n+)2-(A°Asr. 
The splitting between the interacting bands is 

A+_ = 2D = 2Vi2(n- - n+) 



(109) 



(110) 



and the chemical potential is determined implicitly by Eq. (105). 

At zero temperature, these equations can be easily solved. First, suppose n"*" = 0. Then Eqs. 
that 



10^) and (IQC imply 



Sod — 



V12N, 



1/2 



(111) 



Eq. (|11QD becomes A_| = V i2Ns, and the band-filling constraint yields Ng = 27Vo/i. This solution is consistent if (1) 



Vi2A^s/2Ag^g > 1 [Eq^lll)], and (2) fi < , which is equivalent to N0V12 > 1/2. These constraints are precisely 

those obtained in Sec. VB as the boundaries of the normal-state instability and are depicted in Fig. |l^, demonstrating 
that n'^ = over the entire range of the antiferromagnetic phase. Thus, only the lowest, "-" bands are occupied in 
this phase. From Eq. ( p_llD , we also observe that Eod rises from zero continuously along the boundary shared with the 
iVi phase [cf. Fig. |lj], indicating a second-order phase transition at zero temperature. Along the boundary with the 
A^2 phase, however, Eod jumps discontinuously to a finite value, showing that this transition is first-oidei at T — 0. 

This interpretation is confirmed by a calculation of the energy change of the system across the transition, which 
also demonstrates the stability of the antiferromagnetic phase relative to the normal phases. From the relations 



2m* 



and 



and Eqs. dzl), (g|), 



-y 

k 



T 



2 1,2 



^0 Noxe{x) 



2m* 



x] T^OA^oY^(a;), 



and (p^), we can write the total energy [Eq. (pfj)] as 
E 



^ (\ c|2 



(112) 
(113) 

(114) 
(115) 



20 



Inserting the subband occupations, eigenvectors, and chemical potential obtained from the solution of the self-energy 
and band- filling equations at T = 0, we find for iVi phase 



for the N2 phase 



Nl 

8Nr 



(1 + NoVr. 



^SAS 



NoA" 



SAS 



1 - N^V, 



12 



and for the antiferromagnetic (AF) phase 



E 



AF 



A 



■^^ S 

4No 



2NoVr. 



^^SAS 



(116) 



(117) 



(118) 



We see immediately from Eqs. (116) and (Ulq) that, in the range of parameter space in which both iVi and AF 



solutions exist, E < E . In addition, the energy at the Ni-AF phase boundary [cf. Fig. |l2| can be seen to 
be continuous and with continuous first derivatives. The second derivative is discontinuous at the p hase boun dary , 
however, showing that the phase transition is second order. Similarly, it can be shown that Eqs. (117) and ( |ll§| ) 
yield E^^ < E^'^ in the region of parameter space where both solutions exist and that the first derivative of the 
energy difference at the phase boundary is discontinuous. Thus, the antiferromagnetic phase is energetically stable 
with respect to the normal phases whenever the broken-symmetry solution exists, and the phase transition is second 
order from the A'^i phase and first order from the N2 phase at zero temperature in mean-field theory. 

At finite temperature, the self-energy equations in combination with the band-filling constraint must be solved 
numerically. The results of these computations can be used to obtain several quantities which characterize the 
antiferromagnetic phase: the transition temperature and the temperature dependence of the order parameter 
— {pi2 + P21 ) ^'^'i specific heat cy ■ The temperature at which the antiferromagnetic transition occ urs is found 
by linearizing the self-energy equations in the off-diagonal self-energy Sod- This procedure allows Eq. ( |l08|) to be 
written as 



Vi2{n -n"*") 



^SAS 



(119) 



and also leads to the result = Ag^g for Ag^g defined by Eq. (|8^). From these two expressions, we de duce that 



A -I = 2AgAs at Tc- Thus, the critical temperature is d etermined by the simultaneous solution of Eqs. (|119| ) and 

(105) with the constraint A^ = 2AgAs- Using Eq. ( |112| ), we find that the Tc equation may be written after some 

algebra in the reduced variables /3c — Ag^g/Tc, x — Ns/2NoAg^^, and y = iVoVi2 as 



In 



„/3,(a;y+l)/2y 



1 



af3a{xy-l)/2y _ I 



(120) 



The reduced variables x and y are just the axes of the phase diagram in Fig. |l^, in which the contours of constant 
fcsTc/ Ag^g are also shown. We note that the critical temperature can be of the order of Ag^g/Zcs, which may be on 
the order of 10 K for double quantum wells of the type shown in Fig. |[ 

Below the transition temperature, the order parameter becomes finite, and we must solve the full non-linear set of 
equations. To relate the self-energy parameters obtained in this way to the staggered spin density in Eq. (|95|), we first 
note that Eq. (^ij) can be used to show that 



k 



Inserting this result into Eq. (|95| ) and applying Eqs. (|9^) and (| 



[Pl2 + P21] 

At zero temperature, this expression reduces to 



-2(n" 



, we find that the staggered magnetization 

Sod 



P2l\ 



/2A 



SAS 



1/2 



(121) 



(122) 



(123) 



21 



Solving the T > se lf-energy and band filling equations for a variety of interaction strengths and substituting the 



results into Eq. (122) yields the curves in Fig. We see that the staggered magnetization rises rapidly from zero 
below Tc and saturates quickly to its T = value. This behavior is generally expected for an order parameter in 
mean-field theory. 

Another quantity of theoretical and possibly experimental interest is the specific heat. The specific heat is pro- 
portional to the second derivative of the free energy with respect to temperature, so this quantity is discontinuous 
at either a first- or a second-order phase transition. This discontinuity is in principle measurable and would provide 
direct evidence of a thermodynamic phase transition occurring in these systems. Actually observing this discontinuity 
in semiconductor devices of the kind we are considering would be difficult, however, due to the low concentration of 
the relevant electrons. The specific heat is computed from Eq. ( ^ ) for a particular temperature with the eigenenergies 
and chemical potential obtained from the self-consistency equations. Using Eq. ( [76[ ) in Eq. (p9[), the specific heat per 
sample area may be written 



which at low temperature reduces to 



cv = yA^^cciVoT, (125) 



whereApcc is the number of occupied subbands: A"occ = 2 in the A^i ancLj4F phases and A'occ = 4 in the A''2 phase. 
Eq. ( |125| ) is just the usual result obtained from a Sommerfeld expansion.Ej 

Since we are aware of no calculations in the literature regarding the specific heat of a paramagnetic electron gas 
in a quantum well structure, we present in Fig. |l^(a) cv for the non-interacting electron gas to use as a comparison 
for the interacting case. At low temperature, the curves naturally resolve themselves int o two groups according to 



whether one or two subbands are occupied at T = 0. This feature follows directly from Eq. (125) and is emphasized in 
the inset to Fig. |l|(a). Observe that Eq. (^25|) is an inadequate description of cy when T is larger than only a small 



fraction of the subband splitting A ^ar- The precise fraction is density-dependent, as is whether the actual specific 



heat is smaller or larger than Eq. (125) predicts. The reason underlying this behavior is that, unlike the metallic 
case, the Fermi energy and the temperature are often comparable in these quantum well structures, invalidating the 
Sommerfeld expansion. At higher temperatures, the form of the specific heat for both one- and two-subband-occupied 
ground states are similar and have a magnitude at fixed temperature which increases monotonically with density. 
When the temperature become s of the order of the energies of the higher subbands, this two-subband description 



breaks down. Of course, Eq. (124) may still be used to compute cy at these temperatures providing the higher 
subbands are included in this equation and the equation determining the chemical potential. 

With the non-interacting specific heat as a baseline, we can now examine cy in the interacting system in both the 
normal and antiferromagnetic phases. Fig. |l^(b) presents the specific heat as a function of temperature calculated 
for a fixed interaction strength. The different curves show the evolution of cy (T) as the density is increased from the 
single-subband-occupied A'^i phase into the antiferromagnetic phase. The specific heat in the A'^i phase is similar to 
the non-interacting plots in Fig. |l6|(a), exhibiting the low-temperature, linear-in-T behavior expected from Eq. ( |125| ) 
which crosses over at some density-dependent temperature fc^T ^ ^sAS approximately constant value. As one 

enters the antiferromagnetic phase, one sees a discontinuity develop in cy at Tc, which signals the phase transition. 
Unlike superconductivity,EZI there is no universal value for this discontinuity due to the large renormalizing effects of 
the temperature on the band structure parameters. Above Tc, the curves resemble their normal-state counterparts. 



while below Tc, the low-temperature result Eq. (125) seems to hold. 

The latter behavior is not universal, however, as seen from Fig. |l6|(c). This figure shows the evolution of the specific 
heat function at a fixed density and increasing interaction. The system is initially in the two-subband-occupied N2 
phase and becomes antiferromagnetic when A'oVi2 > 1/2 [cf. Fig. |l2|. In the normal state, cy{T) has the shape 
expecte d fro m Fig. Eq(a) with a low-temperature slope characteristic of having two subbands occupied (i.e., A'occ = 4 



in Eq. (125)). Increasing the interaction strength has only minor effects until the antiferromagnetic region of the 



phase diagram is reached, at which point cy develops a discontinuity at Tc. At the same time, the slope of the 



low-temperature specific heat drops by a factor of two in accordance with Eq. (125) and the fact that only the lowest 
interacting bands are occupied in the antiferromagnetic state. Above Tc, the specific heat is qualitatively similar 
to the other normal-state curves. Below Tc, and in contrast to Fig. |l^(b), cy oc T at both low temperatures and 
temperatures near Tc but with different slopes. The slope near Tc is a function of the interaction strength and equals 
the low-temperature value for N0V12 = 1; for N0V12 > 1, the slope near Tc is actually less than the low-temperature 
slope. 
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D. Collective Excitations 



In addition to the ground state and thermodynamic properties of the antiferromagnetic phase, it is also important to 
examine its collective excitations. The first indication of the antiferromagnetic phase transition is the disappearance of 
the intersubband spin-density excitations, and, on general theoretical grounds, one would like to know what replaces 
them in the broken-symmetry phase. Moreover, experimental studies of semiconductor heterostructures by inelastic 
light scattering can measure these excitations, and theory should provide some guidance about the expected signatures 
of the new phase. The latter point is particularly important in light of current searches for this phaseHa 

The basis for studying the collective modes in the antiferromagnetic phase is Eqs. ( ^ ) and supplemented by the 
eigenvalues i^^, eigenvectors e'^, and chemical potential determined as in the preceding Subsection. Because the wave 
functions of the interacting quasiparticles mix different subbands and spins, the bubble ^fijcdil) longer diagonal 
in these indices, and so Eq. (^4|) becomes a 16 x 16 matrix equation in subband and spin space. The remaining 
symmetries in the antiferromagnetic phase do not seem to be amenable to decomposing this matrix equation and 
arriving at an analytic solution, which forces us to adopt a numerical approach. We therefore obtain the interacting 
polarizability llab,cd{(l) by direct numerical inversion of Eq. (|4|) at T = and identify the collective modes from the 
condition 



det 



SacSbd - ^%f{q) {Vfe^cd - Vfd.ce)] = 0. (126) 



Since we are interested primarily in the intersubband spin-density modes, we will focus on the intersubband spin-flip 
polarizability 

ninter(<z) S UatMl) fcT (127) 

with /°^* given by Eq. (|7^). Typical results for the spectral function of this polarizability in the antiferromagnetic 
phase, — ImllintcrC'i'), are shown in Fig. |lj. 

In analyzing these figures, it is useful to keep in mind the following facts about the band structure of the antiferro- 
magnetic phase revealed by the analysis of the preceding Subsection: (1) the interacting band structure consists of two 

sets of two degenerate parabolic subbands separated by an energy , (2) only the lower set of interacting subbands 

are occupied, yielding a single Fermi wave vector kp, and (3) the wave functions corresponding to the interacting 
bands are a mixture of the wave functions of the non-interacting subbands. Facts (1) and (2) indicate that the band 
structure is similar to that in the A^i phase, so we expect to see a region of intersubband particle-hole excitations in 
the spectral function similar to those in Fig. |l^. The kinematics of these excitations require that their spectral weight 
start at huj — A_| a.t q — and spread within the boundaries given by 

— (q - kpf <nuj~A+^+fl<—{q + kpf (128) 
for g > and /i measured from the bottom of the lowest interacting subband. In the normal phases, these inter- 



subband excitations are cleanly separated from the intra-subband excitations by the projection in Eq. (127). In 
the antiferromagnetic state, however, fact (3) indicates that such a separation is impossible, leading to the addi- 
tional low-frequency particle-hole continuum present in Fig. |l^. Kinematics again show that this region is defined by 



Q < uj < 2^(9 + ^-f)^ ~ M for q < 2kp. Although not shown in Fig. 17, at larger q > 2kp, the extent of this region 
of particle-hole excitations is defined by 



{q~kFY~^i<UJ<—{q + kFY-^Ji. (129) 
2m* 2m* 

Within the particle-hole continuum, we find one striking feature at low frequenci es w hich appears in Fig. |l^ as a 
dark, linearly dispersing feature at low frequencies. Examining the solutions of Eq. ( |126| ),EZI we find that this feature 
is in fact a Landau-damped collective mode of the system. The polarization of this collective mode is extracted by 
returning to the full interacting polarizability Ilab,cd{<l) and applying Eqs. (|7j)-(^5|) at the wave vectors and frequencies 
lying on the dispersion curve for this mode. A real-space representation of the resulting spin density displacements is 
given in Fig. |l^(b). To interpret these results, we first note that the antiferromagnetic phase for these calculations has 
its spin density oriented along the x direction. From Fig. |lj(b) , we see that the collective mode corresponds to a wave 
of antiferromagnetic spin displacements normal to this orientation, namely in the y direction, traveling transverse to 
the layering direction of the quantum wells. The spin displacements are in opposite directions in different wells, so 
the net effect is a rotation of the total spin density in the xy plane that preserves the antiferromagnetic correlation 
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of the spin density between the wells. Note that there is another collective mode degenerate with this one which 
corresponds to a rotation of the total spin density in the xz plane; this mode is projected out by our choice of /°^* 
[Eq. @]. 

Outside of the partielc-hole continuum, we can look for the undamped collective modes in the same way as in the 
preceding paragraph.EZI We find a single optical excitation whose dispersion is indicated by the thick black line in 
Fig. 0. The polarization of this mode is obtained from Eqs. (|7i|)-(|75|) and is shown in Fig. |lj(c). From this figure 
and our knowledge of the orientation of the spin density in the antiferromagnetic phase, we conclude that the optical 
mode corresponds to a modulation of the magnitude of the spin density which alters neither its direction in space nor 
the antiferromagnetic correlation between the wells. We note that, unlike the low-frequency mode, the optical mode 
becomes so strongly Landau damped once it enters the particle-hole continuum that it is no longer identifiable. 

The two collective modes that appear in our calculatioss— can be understood from general principles of phase 
transitions involving the breaking of a continuous symmetry .Hca In our case, the continuous symmetry is spin-rotation 
or SU{2) invariance, and the extent to which it is broken is quantified by the staggered spin density N [Eq. |9^)]. The 
collective modes in the broken-symmetry phase correspond to the modulation in space and time of either the direction 
or the magnitude of N and are therefore called phase and amplitude modes, respectively. The polarization of the 
collective excitations discussed above unambiguously identify the two low-frequency. Landau-damped excitations as 
phase modes with orthogonal polarizations and the optical excitation as the amplitude mode. This identification is 
strengthened by the dispersion of these excitations: the phase mode should have an energy which vanishes as q — s- 0, 
since ground states with different orientations of N are degenerate, whereas the amplitude mode should possess an 
excitation gap. These expectations are borne out in Fig. 

Further insight into these collective modes can be gained by following the evolution of the intersubband spin-density 
excitations as we change the parameters of our model and move from the normal to antiferromagnetic phases. We 
start in the single-subband-occupied A^i phase at the point marked Ai in Fig. |2[ Increasing the density at fixed 
interaction, we move towards the antiferromagnetic phase (point A2), causing the three degenerate intersubband 



spin-density excitations to soften as shown in Fig. 13, When the q = excitations vanish, the system enters the 
antiferromagnetic phase and the three intersubband SDEs turn into two degenerate phase modes and an amplitude 
mode as seen in Fig. |l^(a). Increasing the density still further (to point B2), the amplitude mode moves to higher 
frequencies while the phase mode is largely unchanged [Fig. |l7|(b)]. 

This behavior is typical of entry into the antiferromagnetic phase, and can be summarized by a plot of the q = 
intersubband excitation spectrum as a function of the model parameters. Fig. |l^. In Fig. |l8|(a), we see the transition 
just described, in which the system is initially in the Ni phase and the density is increased at fixed interaction strength. 
The collapse of the intersubband spin-density-excitation and the emergence of the amplitude mode is clearly seen. We 
also observe that both the splitting between the interacting subbands and the intersubband charge-density excitations 
are continuous and non-zero across the antiferromagnetic transition, but both have a discontinuity in their first 
derivatives. 

The antiferromagnetic transition has a slightly different character when it proceeds from the two-subband-occupied 
N2 side. Starting in the A'2 phase and increasing the interaction strength at fixed density, we see from Fig. |8|(b) that 
the intersubband spin-density excitation softens as before but the amplitude mode in the antiferromagnetic phase 
appears immediately thereafter at finite frequency. The interacting subband splitting and the intersubband charge- 
density excitation are also discontinuous across this phase boundary. These jumps are a consequence of the first-order 
nature of the T — antiferromagnetic transition from the N2 side. Although finite temperatures will probably restore 
continuity to these curves, the large changes indicated may be experimentally observable and could provide strong 
evidence for the transition. 

Whether the changes in the intersubband excitation spectrum we have discussed are observable in, for example, 
inelastic light scattering experiments depends to a great extent on the spectral weight associated with these features. 
To give an idea of the range of intensities involved, we plot in Fig. a cross section at fixed q of the spectral weight 
of the intersubband polarizability from Fig. 0^b). With increasing frequency, peaks associated with the phase mode, 
the amplitude mode, and the intersubband particle-hole continuum are visible, but the spectral weights associated 
with each peak are vastly different with the response dominated by the low-frequency phase mode. Most inelastic 
light scattering measurements are done at small q and moderate frequencies on this scale, so the phase mode may 
be difficult to observe unless a concerted effort is made to look for it. Indeed, the signature of the antiferromagnetic 
phase in conventional light scattering experiments may simply be the apparent absence of all intensity. We also note 
that the polarization of the scattered light relative to the antiferromagnetic ordering direction may affect the observed 
intensities of these modes. Since the ordering direction is arbitrary, this effect may result in a strong variation in the 
observed light scattering spectra after temperature cycling above Tc or between different samples. 
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VI. DISCUSSION AND SUMMARY 



In this paper, we have studied the magnetic instabiHties of semiconductor quantum wells within the local-density 
approximation to density-functional theory and a self-consistent Hartree-Fock theory. To create a consistent picture 
of the results of these calculations, one must realize that these two formalisms supply complementary information. 
The LDA computations are designed to be quantitatively reliable for the normal-state properties of these quantum 
well structures. The self-consistent Hartree-Fock calculation, on the other hand, is only qualitatively reliable, but it 
is able to describe broken-symmetry phases that cannot be studied within the LDA. In particular, there is little point 
in trying to relate the parameters from the self-consistent Hartree-Fock calculation to the LDA results, because the 
former neglects such real-world effects as the distribution of the donor impurities which the latter includes. Hence, 
the LDA calculations should indicate whether or not the transition occurs and suggest the structures and densities 
at which to look for it, and the self-consistent Hartree-Fock calculations should provide information on qualitative 
feature of the resulting antiferromagnetic phase that can assist experimentalists in identifying it. 

The only qualitative point about which the LDA and self-consistent Hartree-Fock calculations disagree is the 
ordering wave vector of the transition from the two-subband-occupied side of the phase diagr am: LDA yields qc = 
kpi — kF2 while the self-consistent Hartree-Fock calculation gives Qc = 0. As mentioned in Sec. |VB| , the discrepancy 
between the two formalisms can be traced back to the number of subbands included in the calculation and not to 
the form of the interaction, which is taken to be independent of wave vector in both cases. Since one expects a 
calculation including more subbands to be more accurate, it is reasonable to conclude that at least some part of 
the true phase diagram would have Qc ^ 0. That not all of the soft-SDE region of the phase diagram would have 
gc 7^ is demonstrated by that fact that the one-subband-occupied SDKs unambiguously soften at qc = in both 
the LDA and self-consistent Hartree-Fock theories. Since the latter theory focuses on inter-well effects and neglects 
intra-well ones, its predictions regarding inter-well properties such as the long-wave-length intersubband spin-density 
excitations may be qualitatively valid even in the Qc ^ phase. This phase would have a non-trivial spin-density 
modulation transverse to the quantum well layering direction, but whether this modulatimi would be of the form 
of a simple spin-density wave or something akin to an antiferromagnetic Skyrmion latticg£j in zero field cannot be 
determined from the present calculations. Future investigations exploring the (?c 7^ phase, and in particular the 
nature of the crossover between the Qc = and Qc ^ phases, could in principle be performed within a generalization 
of the self-consistent Hartree-Fock formalism discussed in this paper. 

With these caveats in mind, let us summarize the primary results of this work. We have presented a TDLA 
calculation which shows that the intersubband spin-density excitations (SDE) in certain coupled double- and wide 
single-quantum-well structures soften completely in a range of densities around the point where the second subband 
begins to populate (10^" — 10^^ cm^^) and in the absence of an external magnetic field. Based on these calculations, we 
have constructed a phase diagram indicating the structures likely to exhibit this instability. We have also computed 
the excitation spectrum measurable by inelastic light scattering near the instability in order to illustrate how the 
SDE softening would appear in these experiments. Since the TDLDA yields_bpl.h_spin- and charge-density excitation 
spectra which are in very good quantitative agreement with experiment Ji3'E£ft3 the softening of the spin-density 
excitations should be observable in the appropriate range of densities. 

In trying to understand this instability, we have explored the possibility of ferromagnetic transitions in double- 
and single-quantum- well structures by including the spin degree of freedom in a density-functional calculation within 
the LSDA. We find that a ferromagnetic transition occurs in the doublc-quantum-well structures which exhibit the 
SDE softening but that the transition occurs at much lower 10^ cm^'^) densities, implying that the SDE softening 
cannot be associated with ferromagnetism. In square single quantum wells, our computations provide evidence for a 
spin-polarized phase of the electron gas which lies between the Wigner crystal and normal phases. The critical density 
for this transition decreases with increasing well width, demonstrating that exchange-correlation effects are stronger 
in lower dimensions, as expected from a simple Hartree-Fock analysis. 

Having failed to identify the SDE-softcned phase within density-functional theory, we turned to a simple model 
of coupled double quantum wells which we treated within self-consistent Hartree-Fock theory. This model is able to 
reproduce the SDE softening in its normal state, and the polarization of the soft mode indicates that the softening 
signals the onset of antiferromagnetic order in the spin density between the wells. Extending our calculations into 
the antiferromagnetic phase, we find that this phase exists and is stable over a wide range of parameters and that the 
mean-field transition temperature can be of the same order as the symmetric-antisymmetric splitting. In addition, 
we find that the transition to this phase at zero temperature is second-order from the single-subband-occupied side 
of the phase diagram, but first-order from the two-subband-occupied side. Due to the absence of an energy gap 
in the single-particle spectrum, we do not expect strong anomalies in the transport properties to accompany the 
transition; however, our calculations of the electronic specific heat show that, if this quantity is measurable, it will 
show a characteristic discontinuity at the transition temperature. A means of searching for the transition which is 
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more likely to succeed is the measurement of the collective spin-density excitations through inelastic light scattering. 
By computing the spectrum of these excitations in the antiferromagnetic phase, we identify a Landau-damped phase 
mode of the order parameter and a true optical collective excitation corresponding to the amplitude mode. The 
spectral weight associated with the phase mode is large, suggesting that inelastic light scattering experiments should 
look at low frequencies for this characteristic excitation of the antiferromagnetic phase. 

In closing, we note that at least one experimental group has investigated the possibility of an antiferromagnetic 
phase transition of the tyne we describe by performing resonant inelastic light scattering measurements on double- 
quantum-well structures.E3 The results, however, have been mixed. In zero field, the complete softening of the 
spin-density excitation does not seem to appear in the electron density regime (« 5 x 10^" cm~^) predicted by the 
TDLDA theory. It is possible, of course, that impurity-scattering-induced broadening effects make it impossible to 
observe the complete softening of the spin-density excitation. Alternatively, the TDLDA theory may overestimate the 
density range in which the transition occurs, implying that the actual antiferromagnetic instability may take place 
at lower electron densities. A third possibility is that the TDLDA approach may simply be inadequate for studying 
semiconductor quantum wells at the low densities involved. 

More promising are the experimental results in small but finite magnetic fields along the layering direction, which 
do indicate a softening of the intersubband spin-density excitations at the filling factors v = 2 and 6E3 The general 
observation of the softening of the intersubband spin-density excitations in the presence of a magnetic fielcEj is 
qualitatively consistent with the prediction of our zero-field theory, since a magnetic field weak enough so that the 
system in not completely spin-polarized enhances the effects of the interaction by reducing the kinetic energy through 
Landau quantization. Thus, the magnetic field naturally enhances the potential for the type of spin instabilities 
discussed in this paper to appear. Moreover, the basic Hartree-Fock theory underlying the description of the resulting 
broken-symmetry states do not change, although the effects of Landau-level quantization should be included. Such a 
generalization of the self-consistent Hartree-Fock approach is straightforward and is left for future work. 
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FIG. 1. Typical coupled double-quantum-well structure, 
and its self-consistent LDA subband energy levels Ei, 
eigenfunctions (pi, electron density n{z), and Fermi en- 
ergy Ep. Also shown are the effective, Hartree, and ex- 
change-correlation potentials Veff, Vh, and Vxc- The sheet 
density iV^ = 2.68 x 10^^ cm'^ 

FIG. 2. Many-body diagrams used to compute (a) the 
self-energy E, (b) the contribution of the interactions to 
the energy _Eint ~ (-ffint), (c) the generalized polarizability 
n, and (d) the vertex function 7 within the self-consistent 
Hartree-Fock approximation. The solid lines represent 
dressed electronic propag ators [Eq. @] and the dashed lines 
the effective interaction V [Eq. (f43|)], both of which are matri- 
ces in subband and spin space. In order to treat both spin and 
charge polarizabilities with the same equations, the polariz- 
ability is not separated into reducible and irreducible parts. 

FIG. 3. Calculated intersubband charge-density excitation 
-BcDE, spin-density excitation i?sDE, and single-particle exci- 
tation -EsPE = AsAS energies as functions of the 2D electron 
density A^s for a DQW structure with barrier width ds = 
40 Aand well width dw = 139 A. The critical density for the 
instability Ac ~ 0.69 x 10^^ cm~^. The bottom figure shows 
an expanded density range making obvious the re-entrance of 
the normal phase at very low electron density. 



FIG. 7. Energy of the intersubband spin-density excita- 
tions -EsDE as a function of wave vector g in a coupled dou- 
ble-quantum-well system with a barrier width of 40 Aand well 
widths of 139 Afor sheet densities A^s in units of 10^^ cm"'^ 
approaching the critical density Ac — 0.686 x 10^^ cm~^ from 
above (thick lines). The thin lines show the lower boundary 
of the particle-hole continuum, above which the collective ex- 
citations are Landau damped. 

FIG. 8. Calculated Raman scattering spectra in the 
cross-polarization geometry for a double- quantum- well struc- 
ture with a 40 Abarrier width and 139 Awell widths. The 
curves illustrate the signatures of the excitonic instability (a) 
as the sheet density Ns is lowered to the critical density and 
(b) as the wave- vector transfer q is increased at constant sheet 
density A^s = 0.7 x lO" cm"^ 

FIG. 9. Typical wide square quantum well given by the 
bare confining potential Vconf, and its self-consistent LDA 
subband energy levels E„, eigenfunctions </!>„, electron den- 
sity n{z), Fermi energy Ep, and efi'ective, Hartree, and ex- 
change-correlation potentials Veff, Vh, and Vxc- The sheet 
density is A^s = 0.9 x 10^^ cm~^. The figure shows how the 
electronic density profile becomes localized on the sides of the 
well, similar to the profile in a double quantum-well system. 
Bottom: lowest energies in expanded scale. 



FIG. 4. Dependence on the sheet density As of (a) 
the mean-field vertex correction | (7x012X12 I (solid line) 
and the absolute value of the lowest-order polarizabil- 
ity I X12 I (dashed line) and (b) the spin-polarized ex- 
change-correlation-induced vertex correction f/xci2 for the 
double-quantum-well structure in Fig. |l| Note that the elec- 
tron gas is unstable the range of densities where | UxcXi2 |> 1 
(see text). 

FIG. 5. Calculated zero-temperature phase diagram for 
double quantum wells in terms of the sheet density As and 
the barrier width ds for fixed well widths dw = 139 A. Inset: 
phase diagram for fixed d_g = 40 Ain terms of the sheet den- 
sity and well widths. The normal (N) and the triplet excitonic 
(E) phases are shown. 



FIG. 10. Calculated intersubband charge-density excita- 
tion energy Ecde, spin-density excitation energy Esde, and 
single-particle excitation energy Espe = As as as func- 
tions of the 2D electron density As for the wide square- well 
structure shown in Fig. ^ Note the collapse of -Esde for 
As ~ 0.2 - 0.4 X 10" cm~^ 

FIG. 11. Approximate zero-temperature, spin-polarization 
phase diagram of single square wells calculated in the lo- 
cal-spin-density approximation in terms of the sheet density 
As and well width ds (top) and the Ts parameter and well 
width (bottom). 



FIG. 6. Calculated zero-temperature phase diagram for 
coupled double quantum wells in 

terms of = (ttNs)~^^^ /ao and the dimensionless symmet- 
ric-antisymmetric subband splitting Asas /{e^ /^{dw + ds)), 
where ao and e are the Bohr radius and dielectric constant for 
GaAs and ds {dw) is the barrier (well) width. Solid circles 
correspond to dw = 139 Aand various ds and are taken from 
the main part of Fig. |^, while the crosses correspond to ds 
— 40 Aand various dw and are obtained from the inset to 
Fig. |5[ The normal (N) and the triplet excitonic (E) phases 
are shown. 
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FIG. 12. Mean-field phase diagram of the antiferromag- 
netic sector of the point-contact model described in the 
text. The independent variables are the normalized intersub- 
band interaction matrix element N0V12 and the sheet density 
Ns /2NoA^j^g, where A'o ~ m* in the single-spin den- 
sity of states and AgA3 is the splitting between the lowest 
two subbands when V12 = 0. The other interaction matrix 
elements Vu = V22 ~ 0. The regions correspond to the nor- 
mal (paramagnetic) phase with one subband occupied (A''i), 
the normal phase with both subbands occupied {N2), and the 
antiferromagnetic phase (AF). Contours in the antiferromag- 
netic region of the phase diagram are the computed values of 
the critical temperature Tc for the antiferromagnetic transi- 
tion in units of Ag^g/fc^. Observe that ksTc can be larger 
than Ag^g, indicating that the antiferromagnetic phase may 
persist to observable temperatures. The other labels in the 
figure identify points for future reference. 



FIG. 13. Disper- 

sion of the intersubband spin-density excitations (SDEs) as 
the sheet density Ns approaches the antiferromagnetic phase 
from the low-density, one-subband-occupied side computed in 
the antiferromagnetic sector of the point-contact interaction 
model discussed in the text with N0V12 = 1.0. The thick 
lines show the energy hu) of the SDEs in units of the renor- 
malized splitting of the lowest two subbands Ag^s as a func- 
tion of wave vector q relative to the Fermi wave vector kp- 
The thin lines show the boundaries of the particle-hole con- 
tinuum, within which the collective excitations are damped. 
The sheet densities and corresponding points in the phase 
diagram of Fig. ^ are given in the figure. The transition 
to the antiferromagnetic phase occurs when the g = SDEs 
soften at A^s/2A'oAgAs ~ 1.0. A similar softening appears 
in Fig. 1^, which shows the approach to the antiferromagnetic 
phase from the high-density side computed within the LDA. 



FIG. 14. Polarization of (a) an intersubband spin-density 
excitation in the normal phase and of the (b) phase 
(Nambu-Goldstone) and (c) amplitude modes in the antifer- 
romagnetic phase with the spin density oriented along the x 
direction computed as described in the text. The configura- 
tion of the quantum wells is as in Fig. |^. The two planes are 
sections through this geometry normal to the layering direc- 
tion and are located in the center of each well. The distances 
in these planes are measured in units of the wavelength of 
the collective excitation A, whose propagation is in the x di- 
rection. The arrows show the direction and magnitude of the 
spin density modulation induced by the collective excitations. 
These modes have the form of a traveling wave, so the spin 
modulation at a different time is obtained by shifting these 
pictures along the x direction. Since the total spin density is 
the sum of the antiferromagnetic polarization and the modu- 
lations shown in (b) and (c) , the identification of these modes 
with the phase and amplitude motions is apparent. 



FIG. 15. Staggered spin density {pi2 + P21), which is the 
order parameter for the antiferromagnetic phase transition 
discussed in the text [cf. Eq. (^6|)], normalized by the elec- 
tron sheet density Ns as a function of temperature T in 
units of the splitting between the two lowest subbands in 
the non-interacting limit Ag^s- The curves are computed for 
A''a/2iVoAsAs = 2.5 and iVoV'12 = 0.55 to 0.80 in increments 
of 0.05; the lower and upper values correspond to points Di 
and D2 in the phase diagram of Fig. ^ Inset: Expectation 
value of the spin density (5^(2;)) in real space as a function 
of the distance along the layering direction z for the double 
quantum well of Fig. ^ in the antiferromagnetic phase. Note 
that (3^(2)) = in the paramagnetic phase. 



FIG. 16. Electronic specific heat at constant volume cv 
in normalized by Nq/^^^c^/Z as a function of the tem- 
perature T in units of AgAs/fcs for (a) A'oVi2 = and 
Ns/2No/^Ia_s = 0.1 to 2.0 in increments of 0.1, (b) iVoVi2 = 1 
and A^s/2A'oAgAg ~ 0.5 to 1.5 in increments of 0.1 (i.e., along 
the line from to B2 in Fig. |l|), and (c) iVs/2iVoAgAS = 
2.5 and A*'oVi2 = 0.2 to 0.8 in increments of 0.05 (i.e., along 
the line from Ci to D2 in Fig. ^). The inset in (a) enlarges 
the low-temperature portion of the main figure in order tosee 
the deviation from the analytic low-T expression, Eq. (125). 
Entrance into the antiferromagnetic phase is signaled by a 
discontinuity in the specific heat that is apparent in (b) and 
(c) but absent in (a). 



FIG. 17. Dispersion of the intersubband spin density col- 
lective modes in the antiferromagnetic phase at the points 
(a) B\ and (b) B2 of the phase diagram of Fig. FL3, corre- 
sponding to (iV,/2AroAgAs, A^oVi2) = (1.1,1.0) and (1.5,1.0), 
respectively. The thick lines are the energy Ttw in units of 
the non-interacting subband splitting Ag^g of the amplitude 
mode of the antiferromagnetic order parameter as a function 
of the wave vector q in units of q\ — m* Ag^g/ft^. The shaded 
region is the particle-hole continuum, with darker shades rep- 
resenting larger spectral weight than lighter shades on a log- 
arithmic intensity scale. The dark linear feature is the phase 
or Nambu-Goldstone mode of the order parameter. Note that 
this mode is damped by particle-hole excitations and that in- 
tra-subband excitations enter into the spectrum due to the 
mixing of the non-interacting wave functions in the symme- 
try-broken phase. 
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FIG. 18. Normalized energy ftc^/A^As of the g = 0, T = 
collective excitations and interacting subband splitting for (a) 
a fixed interaction N0V12 = 1.0 and varying sheet density 
Ns /2NoAgj^g and (b) varying interaction A'oVi2 and fixed 
sheet density A^s/2A'oAgAs = 2.5 in the antiferromagnetic 
sector of the point-contact model discussed in the text. Il- 
lustrated are the intersubband spin-density (SDE, solid lines) 
and charge density excitations (CDE, dotted lines), the renor- 
malized subband splitting Ag^s (SPE, dot-dashed lines), and 
the amplitude mode in the antiferromagnetic phase (dashed 
line). The top axes show the corresponding points in the 
phase diagram in Fig. As seen in (b), the collective mode 
energies are discontinuous across the two-subband-occupied 
{N2) to antiferromagnetic phase boundary at T = 0, indicat- 
ing a first-order transition. 



FIG. 19. Imaginary part o f the intersubband spin-flip po- 



larizability — Imllintcr [Eq. (127)] as a function of the ex 



citation energy hw relative to the non-interacting subband 
splitting AgAs for q — OAqa at the point B2 in the an- 
tiferromagnetic region of the phase diagram of Fig. |l^ [see 
also Fig. |l7|(b)]. This quantity is related to the intensity of 
the intersubband response in inelastic light scattering exper- 
iments. The phase or Nambu-Goldstone mode dominates the 
low-frequency response, the amplitude mode shows much less 
intensity, and the particle-hole excitations have a very weak 
signal. This spectrum is computed using a finite scattering 
rate 7 — COlAg^s to simulate the effects of impurities. 
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